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ABSTRACT 

This  report  is  the  written  version  of  a series  of  lectures 
presented  by  the  author  at  a NATO  advanced  Studies  Institute 
held  at  the  University  of  Liege,  Belgium,  during  the  period 
September  6-17,  1976.  The  report  presents  an  exposition  of  the 
results  of  the  author's  research  on  acoustic  wave  propagation 
during  the  last  five  years.  The  principal  result  of  this  research 
is  the  "method  of  spectral  and  asymptotic  analysis"  for  determining 
the  structure  of  acoustic  waves  in  unbounded  media.  In  this 
report  the  method  is  explained  and  illustrated  by  its  application 
to  the  following  four  classes  of  acoustic  wave  propagation  problems. 

1.  Scattering  by  bounded  obstacles  in  a homogeneous  unlimited 
fluid. 

2.  Propagation  and  scattering  in  simple  and  compound  tubular 
waveguides. 

3.  Propagation  in  plane  stratified  fluids. 

4.  Propagation  in  crystals. 
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SPECTRAL  AND  ASYMPTOTIC  ANALYSIS  OF  ACOUSTIC  WAVE  PROPAGATION 


Calvin  H.  Wilcox 

Department  of  Mathematics,  University  of  Utah, 
Salt  Lake  City,  Utah,  USA 


1.  INTRODUCTION 

Classical  theories  of  acoustic  wave  propagation  provide  a 
wealth  of  examples  of  boundary  value  problems  for  evolution  partial 
differential  equations.  These  problems  may  be  described  categori- 
cally as  initial-boundary  value  problems  for  certain  systems  of 
linear  hyperbolic  partial  differential  equations  with  variable 
coefficients.  However,  the  known  existence,  uniqueness  and  regu- 
larity theorems  for  these  problems  are  only  a first  step  toward 
understanding  the  structure  of  the  solutions.  To  obtain  a deeper 
insight  it  is  essential  to  discover  how  the  nature  of  the  solutions 
changes  with  the  geometry  of  the  boundary  and  with  the.  coefficients. 
An  examination  of  recent  scientific  literature  on  acoustics  reveals 
a great  variety  of  physically  distinct  phenomena.  Examples  include 
phenomena  associated  with  acoustic  wave  propagation  in  stratified 
fluids,  anisotropic  solids  such  as  crystals  and  man-made  composites, 
open  and  closed  waveguides,  periodic  media  and  many  others.  A 
theory  which  treats  all  of  these  phenomena  on  the  same  footing  can 
provide  only  the  most  superficial  information  about  the  structure 
of  acoustic  waves. 

The  purpose  of  these  lectures  is  to  present  a method  for 
determining  the  structure  of  acoustic  waves  in  unbounded  media. 

The  method  will  be  explained  in  the  context  of  four  specific  classes 
of  propagation  problems.  No  attempt  will  be  made  to  formulate  the 
most  general  problem  that  can  be  analyzed  by  the  method.  Indeed, 
such  a formulation  wouLd  necessarily  be  too  abstract  to  be  useful. 
However,  it  will  be  clear  from  the  examples  that  the  method  is 
applicable  to  many  other  wave  propagation  problems,  both  in  acous- 
tics and  in  other  areas  of  physics. 
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It  will  be  helpful  to  outline  the  main  steps  of  the  method 
here  before  passing  to  a detailed  discussion  of  specific  cases. 

The  method  is  based  on  the  fact  chat  the  states  of  an  acoustic 
medium  which  occupies  a spatial  domain  ft  C R3  can  be  described  by 
the  elements  of  a Hilbert  space  JC  of  functions  on  ft.  The  evolution 
of  an  acoustic  wave  in  the  medium  is  then  described  by  a curve 
t u(t,»)  € j(.  Moreover,  there  is  a selfadjolnt  real  positive 
operator  A on  Jf,  determined  by  the  geometry  of  ft  and  the  physical 
properties  of  the  medium,  such  that  the  evolution  of  acoustic  waves 
in  the  medium  is  governed  by  the  equation 

2 

|~  + Au  = 0 (1.1) 


It  follows  that  the  evolution  is  given  by 

u(t,*)  - Re  {exp  (-itA1^2  )h)  (1.2) 

where  h £ 3f  characterizes  the  initial  state  of  the  wave. 


The  spectral  theorem  may  be  used  to  construct  the  solution 
operator  exp  (-itA1^2) . However,  the  very  generality  of  this 
theorem  implies  that  it  can  give  little  specific  information  about 
the  structure  of  the  wave  functions  u(t,x).  Accordingly,  the  next 
step  in  the  method  is  to  construct  an  eigenfunction  expansion  for 
A.  In  each  of  the  cases  discussed  below  A has  a purely  continuous 
spectrum  and  the  eigenfunctions  are  therefore  generalized  eigen- 
functions. They  define  a complete  set  of  steady-state  modes  of 
propagation  of  the  medium  and  the  most  general  time-dependent 
acoustic  wave  in  3C  can  be  constructed  as  a spectral  integral  over 
these  modes. 


The  final  step  in  the  method  is  an  asymptotic  analysis  for 
t 00  of  the  spectral  integral  representing  u(t,x).  The  result  is 
an  asymptotic  wave  function  u°°(t,x)  which  approximates  u(t,x)  in 
3f  when  t -*•  that  is. 


lim  fl  u(t,  •) 

£-K)0 


- 0 


(1.3) 


Stronger  forms  of  convergence  can  also  be  proved  under  appropriate 
supplementary  hypotheses  about  the  medium  and  its  Initial  state. 


The  result  (1.3,)  offers  a fundamental  insight  into  the  nature 
of  transient  acoustic  waves  in  unbounded  media.  For  it  is  found  in 
each  case  that  the  form  of  the  asymptotic  wave  function  u (t,x)  is 
determined  entirely  by  the  geometry  of  the  domain  ft  and  the  physical 
characteristics  of  the  medium  that  fills  it.  Only  the  fine  struc- 
ture of  u“'(t,x)  depends  on  the  initial  state  of  the  wave.  Thus  in 
the  simple  case  of  a homogeneous  fluid  filling  R3,  u°°(t,x)  is  a 
spherical  wave: 


u°°(t,x)  = F (r  - t,0)/r,  x = r0,  1 0 [ = 1 
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(1.4) 


The  initial  state  affects  only  the  shape  of  the  profile  F(x,0).  In 
other  cases  the  form  of  u (t,x)  is  entirely  different,  but  in  each 
case  the  form  of  u (t,x)  is  determined  solely  by  the  geometry  and 
physical  characteristics  of  the  medium.  In  each  case  u“(t,x)  gives 
the  final  form  of  any  transient  wave  in  the  medium.  The  details  of 
how  the  wave  is  excited  have  only  a secondary  effect  on  the  ultimate 
waveform . 

The  remainder  of  these  lectures  is  organized  as  follows.  The 
fundamental  boundary  value  problems  of  acoustics  are  formulated  in 
section  2.  The  spectral  and  asymptotic  analysis  of  the  four  classes 
of  propagation  problems  is  presented  in  sections  3 through  8.  The 
four  classes,  which  are  physically  quite  different,  were  chosen  to 
illustrate  the  flexibility  and  scope  of  the  method.  In  each  of  the 
four  classes  there  is  a special  case  for  which,  because  of  addi- 
tional symmetry,  the  eigenfunctions  can  be  constructed  explicitly. 
The  remaining  cases  of  the  class  are  then  treated  as  perturbations 
of  the  special  case.  When  used  in  this  context,  perturbation  theory 
is  usually  called  the  steady-state,  or  time-dependent,  theory  of 
scattering.  The  first  class  of  problems  treated  below  corresponds 
physically  to  the  scattering  of  acoustic  waves  by  bounded  obstacles 
immersed  in  a homogeneous  fluid.  Mathematically,  it  is  an  initial- 
boundary problem  for  the  d'Alembert  equation  in  an  exterior  domain 
ft  C R3  (R3  -ft  compact) . The  simple  special  case  where  ft  = R3  is 
treated  in  section  3 and  the  general  case  in  section  4.  The  second 
class  of  problems  deals  with  tubular  waveguides.  Thus  ft  is  the 
union  of  a bounded  domain  and  a finite  number  of  semi-infinite  cyl- 
inders. The  special  case  of  a single  cylinder  is  treated  in  section 
5 and  the  general  case  in  section  6.  The  third  class  of  problems, 
treated  in  section  7,  deals  with  acoustic  wave  propagation  in  plane 
stratified  fluids  filling  a half-space.  Here  the  novel  feature  is 
the  possibility  of  the  trapping  of  waves  by  total  internal  reflec- 
tion. The  fourth  and  final  class  of  problems,  dealing  with  acoustic 
waves  in  crystalline  solids,  is  discussed  in  section  8.  The  new 
feature  in  this  case  is  the  anisotropy  which  has  a profound  effect 
on  the  form  of  the  asymptotic  wave  functions. 

The  results  presented  below  are  based  primarily  on  the  author's 
research.  Sections  3 and  4 are  based  on  the  author's  monograph  on 
"Scattering  Theory  for  the  d'Alembert  Equation  in  Exterior  Domains" 
[42].  The  spectral  theory  of  acoustic  wave  propagation  and  scatter- 
ing in  tubular  waveguides  was  developed  by  C.  Goldstein  [9-12]  and 
by  W.  C.  Lyford  [21,22].  More  recently,  J.  C.  Guillot  and  the 
author  [13]  have  developed  the  theory  for  domains  ft  which  are  the 
union  of  a bounded  domain  and  a finite  number  of  cylinders  and 
cones.  Sections  5 and  6 present  spectral  and  scattering  theory  for 
tubular  waveguides  following  the  plan  of  [13]  . Sections  7 and  8 are 
based  on  the  author's  publications  [39,40,43,44]. 
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The  goal  of  these  lectures  Is  to  provide  an  introduction  to 
the  method  of  spectral  and  asymptotic  analysis  of  wave  propagation. 
Therefore,  the  lectures  emphasize  concepts  and  results,  rather  than 
techniques  of  proof.  Proofs  of  the  results  given  here  may  be  found 
in  the  references  listed  at  the  end  of  the  lectures. 


2.  BOUNDARY  VALUE  PROBLEMS  OF  ACOUSTICS 

Acoustic  waves  are  the  mechanical  vibrations  of  small  amplitude 
that  are  observed  in  all  forms  of  matter.  The  classical  equations 
of  acoustics  are  the  linear  partial  differential  equations  which 
govern  small  perturbations  of  the  equilibrium  states  of  matter. 
Derivations  of  these  equations  from  the  laws  of  mechanics,  together 
with  a discussion  of  their  range  of  validity,  may  be  found  in 
[3,4,8,20,31].  In  this  section  the  equations  and  their  physical 
interpretation  are  reviewed  briefly  and  the  principal  boundary 
va'  ue  problems  for  them  are  formulated  and  discussed.  Applications 
of  the  equations  to  particular  classes  of  acoustic  wave  propagation 
problems  are  developed  in  sections  3 through  8. 

The  following  notation  is  used  throughout  the  remainder  of  the 
lectures.  t (E  R denotes  a time  coordinate,  x - (x1,x2,x3)  € R3 
denote  Cartesian  coordinates  of  a point  in  Euclidean  space.  ft  c R3 
denotes  a domain  in  R3  and  3ft  denotes  the  boundary  of  ft. 
v = Oj.Vj.Vj)  = v(x)  denotes  the  unit  exterior  normal  vector  to  3ft 
at  points  x £ 3ft  where  it  exists.  The  equations  of  acoustics  are 
written  below  in  the  notation  of  Cartesian  tensor  analysis.  In 
particular,  the  summation  convention  is  used.  Acoustic  waves  in 
fluids  (gases  and  liquids)  and  solids  are  discussed  separately. 

The  simpler  case  of  fluids  is  treated  first. 


2.1  Acoustic  waves  in  fluids 

The  case  of  an  inhomogeneous  fluid  occupying  a domain  ft  C R3 
is  considered.  The  propagation  of  acoustic  waves  in  such  a fluid 
is  gove’iod  by  two  functions  of  x 6 ft: 

p = p(x) , the  equilibrium  density  of  the  fluid  (2.1) 

and 

c = c(x),  the  local  speed  of  sound  in  the  fluid  (2.2) 

The  state  of  the  acoustic  field  in  the  fluid  is  determined  by 

v.  *=  v.(t,x),  the  velocity  field  of  the  fluid  at  (2.3) 

J ' time  t and  position  x 
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and 

p = p(t,x),  the  pressure  field  of  the  fluid  at  (2.4) 

time  t and  position  x 

Moreover,  it  is  assumed  that 

p(t,x)  = p q (x)  + u ( t , x)  (2.5) 

where  p0(x)  is  the  equilibrium  pressure  of  the  fluid  and  u(t,x) 
remains  small.  With  this  notation  the  equations  satisfied  by  the 
acoustic  field  in  the  fluid  are 


— i + — - =0  i = 1 2 3 

3t  p(x)  3Xj 

+ c2(x)  p(x)  = 0 
1 

Elimination  of  the  velocity  field  gives  the  single  equation 

32u  2 / \ / \ 8 1 3u  1 ~ 

at2  " C (x)  p(x)  3x.  p(x)  3x.J  = ° 

3 v 32 


(2.6) 


(2.7) 


(2.8) 


for  the  pressure  increment  u = p - p0 . Moreover,  if  u is  known 
then  the  velocity  field  v^  can  be  calculated  from  (2.6). 

The  wave  equation  (2.8)  roust  be  supplemented  by  a boundary 
condition  at  the  fluid  boundary  30.  Two  physically  distinct  cases 
are  considered  here.  The  first  case  is  that  of  a free  boundary 
30.  Here  the  pressure  at  the  boundary  is  unperturbed;  that  is, 


u ^ = 0 if  30  is  a free  boundary 


( 2 - 9 ) 


This  condition  is  often  used  to  represent  an  air-water  interface 
in  the  theory  of  underwater  sound.  The  second  case  is  that  of  a 
rigid  boundary  30.  Here  the  normal  component  of  the  fluid  velocity 
must  vanish:  v.v.  = 0 on  30.  It  follows  from  (2.6)  that 


= v.  -5-—  = 0 if  30  is  a rigid  boundary 

J 3xj  30 


(2.10) 


The  solvability  of  the  boundary  value  problems  (2.8),  (2.9)  and 
(2.8),  (2.10)  is  discussed  below,  after  the  discussion  of  acoustic 
waves  in  solids. 
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2.2  Acoustic  waves  in  solids 


The  case  of  an  inhomogeneous  elastic  solid  occupying  a domain 
S2  C R3  is  considered.  The  propagation  of  acoustic  waves  in  such  a 
solid  is  governed  by  the  following  functions  of  x £ 11: 

p = p(x),  the  equilibrium  density  of  the  solid  (2.11) 

and 

c„k  = c^(x),  the  stress-strain  tensor  for  the  solid  (2.12) 

Jim  Jim 

The  stress-strain  tensor  must  have  the  symmetry  properties  [4]. 


clm  = c£m  = cm£  = ckj  for  a11  j’k’£’m  = X’2’3  (2>13) 

It  follows  that  the  81  components  4m(x>  are  determined  by  21 
functions.  The  state  of  the  acoustic  field  in  the  solid  is 
determined  by 


u 


j 


u.(t,x),  the  displacement  field  of  the  solid 
3 at  time  t and  position  x 


(2.14) 


and 


0.,  = O , (t,x),  the  stress  tensor  field  of  the 


jk  ~jk 


solid  at  time  t and  position  x 


(2.15) 


Moreover,  the  stress  tensor  field  is  symmetric: 

a..  = a,  . for  all  j,k  = 1,2,3  (2.16) 

jk  kj 

With  this  notation  the  equations  satisfied  by  the  acoustic  field 
in  the  solid  are 


V 

Jim,  , . , . , , 

= cik(x)  ir'  • J’k  " 1’2’3 

(2.17) 

in 

L = ^ i - 1 2 3 

p(x)  3xk  ’ 3 

(2.18) 

3t2 

Elimination  of  the  stress  tensor  gives  the  equations 


!luj.  _ 

3t2  p(x)  3x. 
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Jim.  . 

cjk(x)  3T 

J m 


- 0,  j = 1,2,3 


(2.19) 


/ 


1 


for  the  displacement  field  Uj . Moreover,  if  uj  is  known  then 

the  stress  tensor  field  a..  can  be  calculated  from  (2.17). 

Jk 

The  wave  equation  (2.19)  must  be  supplemented  by  boundary 
conditions  at  the  boundary  DO  of  the  solid.  Only  the  cases  of 
free  and  rigid  boundaries  will  be  considered  here.  In  the.  first 
case  the  normal  component  of  the  stress  must  vanish  at  the 
boundary.  Hence 


°3k\ 


DO  cjk  Dx  Vk 

m 


= 0 if  DO  is  a free 
DO  boundary 


(2.20) 


In  the  second  case  the  displacement  must  vanish  at  the  boundary; 
that  is. 


u . 
J 


DO 


= 0 if  DO  is  a rigid  boundary 


(2.21) 


2.3  Energy  integrals 


One  of  the  most  important  formal  properties  of  the  equations 
of  acoustics  is  the  existence  of  quadratic  energy  integrals.  The 
first  order  system  (2.6),  (2.7)  for  acoustic  waves  in  fluids  has 
the  quadratic  energy  density 


n(t,x)  = y< 


pMv.vj  + -!(-)S(kT 


u2  t> 


and  corresponding  energy  integral 


(2.22) 

(2.23) 

where  dx  = dxjdxjdxj  denotes  Lebesgue  measure  in  R3 . The  energy 
density  for  the  derived  field  v!  = Dvj/Dt,  u'  = Du/Dt,  which  also 
satisfies  the  field  equations  (2.6),  (2.7),  can  be  written 


E(v1,v2,v3,u,K,t) 


n(t,x)dx 


n'(t,x)  = f < 


1 Du  Du  1 [Pul 

p(x)  Dx.  Dx.  c'  (x)p(x)  [DLj 


2) 


by  (2.6).  The  integral 


E(u,K, t) 


n'  (t,x)dx 


(2.24) 


(2.25) 


K 


is  an  energy  integral  for  solutions  of  the  scalar  wave  equation 
(2.8).  The  importance  of  these  integrals  in  the  theory  of  acoustic 
waves  derives  from  the  conservation  laws  for  them.  In  differential 
form  they  state  that 
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(2.26) 


and 


Mt.i.xl  = - — - (uv.) 
9t  9x.  3 


9n* (t,x)  _ 9 

9t  9x . 

3 


9u  9u 


P (x)  9t  9xj  J 


(2.27) 


These  equations  follow  immediately  from  (2.6),  (2.7)  and  the 
definitions.  The  integral  forms  of  the  conservation  laws  follow 
from  (2.26),  (2.27)  and  the  divergence  theorem.  They  may  be 
written 


dE(v3 »v2,v3,u,K,t)/dt  = ~ 


and 


dE(u,K, t) / dt  = 


9K 

1 9u  9u 


u(v.v .)dS 
3 y 


9K 


p(x)  9t  9v 


dS 


(2.28) 


(2.29) 


where  K C R3  is  a..y  domain  for  which  the  divergence  theorem  is 
valid  and  dS  is  the  element  of  area  on  9K.  In  particular,  if 
u(t,x)  is  a solution  of  (2.8)  which  satisfies  (2.9)  or  (2.10)  then 
(2.29)  implies  that  dE(u,^,t)/dt  = 0. 

The  equations  for  acoustic  waves  in  solids  have  an  analogous 
quadratic  energy  integral 


E(u1,u2,u3,K,t)  = 
with  density 


k 


n(t,x)dx 


n(t,x)  = ~ < 


p(x) 


9ui  9u3  £m . . 9u£ 

9t  9t  + Cjl<-(X  3^  9xm 


(2.30) 


(2.31) 


The  corresponding  conservation  law,  which  follows  from  (2.19),  is 

o 3u„  3u. 

l2-32> 
m 

in  differential  form  and 


9n(t,x)  = 9 

9t  9x^ 


dE(Uj ,u2 ,u3 ,K, t) /dt  = 


9K 


9u 

(0jkV  *9t  ds 


(2.33) 


in  integral  form.  In  particular,  solutions  of  (2.19)  which  satisfy 
(2.20)  or  (2.21)  also  satisfy  dE (u3 , u2 , u3 , t) /dt  = 0. 
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The  preceding  remarks  emphasize  the  mathematical  relationship 
of  the  quadratic  energy  integrals  to  the  field  equations  of  acous- 
tics. The  term  "energy"  has  been  used  because  in  certain  cases 
the  integrals  can  be  interpreted  as  the  portion  of  the  energy  of 
the  acoustic  field  that  is  in  the  set  K at  time  t.  This  interpre- 
tation is  not  always  correct  because  the  linear  equations  of 
acoustics  are  only  a first-order  approximation  to  more  complicated 
nonlinear  equations  and  the  energy  densities  defined  above  are 
second-order  quantities.  Hence,  it  is  possible  that  other  second- 
order  terms  which  were  dropped  in  the  linearization  should  be 
included  in  the  energy  densities.  A correct  calculation  of  the 
energy  must  begin  with  the  original  nonlinear  problem.  A discussion 
of  these  problems  may  be  found  in  [8,31]  for  the  case  of  fluids  and 
in  [4]  for  the  case  of  solids. 

It  is  important  to  realize  that  the  energy  integrals  defined 
above  play  an  essential  role  in  the  theory  of  acoustic  fields, 
whether  or  not  they  represent  the  actual  physical  energy  of  the 
fields.  Indeed,  it  was  shown  in  [33]  and  [34]  that  the  existence 
of  these  integrals  implies  the  existence  and  uniqueness  of  solu- 
tions to  the  basic  initial-boundary  value  problems  for  acoustic 
fields.  Moreover,  recent  work  on  eigenfunction  expansions  and 
scattering  theory  makes  use  of  Hilbert  spaces  based  on  energy 
integrals.  The  one  indispensible  hypothesis  that  must  be  made  is 
that  the  quadratic  forms  (2.22)  or  (2.24)  and  (2.31)  be  positive 
definite.  For  (2.22)  and  (2.24)  this  means  that 

p(x)  > 0 and  c2(x)  > 0 for  all  x € ft  (2.34) 

In  any  case,  these  hypotheses  are  essential  because  of  the  physical 
interpretation  of  p(x)  and  c(x).  The  form  (2.31)  is  positive 
definite  if  p(x)  > 0 and 

cjk(x)  C£m  Cjk  > 0 for  a11  x G a and  = f 0 (2'33) 

The  last  condition  can  also  be  expressed  by  means  of  the  well-known 
determinantal  criteria  for  a quadratic  form  to  be  positive  definite. 
It  is  assumed  throughout  these  lectures  that  (2.34)  and  (2.35)  are 
satisfied. 

It  has  been  shown  that  the  acoustic  fields  in  both  fluids  and 
solids  satisfy  partial  dif ferential  equations  of  the  form 


+ Au  = 0 (2.36) 

where  A is  a second  order  partial  differential  operator  in  the 
space  variables.  In  the  case  of  fluids  u(t,x)  S R,  Au(t,x)  € R 
and 
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Au 


-c'(x) 


1 9u  ] 

p(x)  9x^  J 


(2.37) 


while  in  the  case  of  solids  u(t,x)  = (u, (t,x) ,n2(t,x) ,u3(t,x))  € R3 , 
Au(t,x)  e R3  and 


(Au) 


1 9_ 

p(x)  9xk 


£m 

c 

Jk 


(x) 


j = 1,2,3 


(2.38) 


Thus  in  both  cases  the  evolution  of  acoustic  waves  in  a medium 
which  fills  a domain  ft  C R3  is  described  by  the  solution  of  an 
initial-boundary  value  problem  of  the  form 


3fu 

at2 


+ Au  = 0 for  t > 0,  x € ft 


(2.39) 


8u  = 0 for  t > 0,  x e 9ft  (2. AO) 

u(0,x)  = f(x)  and  9u(0,x)/9t  = g(x)  for  x £ ft  (2.41) 


Here  (2.40)  represents  one  of  the  boundary  conditions  (2.9), 
(2.10)  in  the  case  of  a fluid  and  (2.20),  (2.21)  in  the  case  of  a 
solid. 


It  is  interesting  to  note  that  the  positive  definiteness  of 
the  energy  densities,  hypothesized  above  on  physical  grounds,  im- 
plies the  hyporbolicity  of  the  equation  (2.36).  It  follows  that 
the  initial-boundary  value  problem  (2.39)  - (2.41)  has  compact 
domains  of  dependence  and  influence  [6,33].  In  physical  terms 
this  means  that  acoustic  waves  propagate  into  undisturbed  portions 
of  a medium  with  finite  speed. 


A simple  and  rigorous  solution  theory  for  the  initial-boundary 
value  problem  (2.39)  - (2.41)  can  be  based  on  the  theory  of  self- 
adjoint  operators  in  Hilbert  space.  This  possibility  follows  from 
the  divergence  theorem  which  implies  the  formal  selfadjointness  of 
the  operators  A relative  to  suitable  inner  products.  Indeed,  for 
the  operator  (2.37)  the  divergence  theorem  implies 


Au  v c 2(x)  p_1(x)dx 
Jft 


(2.42) 


9u  9v 


Jn  9x.  9x. 

j j 


p_1 (x)dx  - 


* 9ft 


(x)dS 


and  hence 


(Au  v - u Av}c~2(x)  p_1(x)dx  = 

Jft 


9v 

9v 


(2.43) 
I"  v}  p_1 (x)dS 
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Thus  ir  an 
(u,v) 

then 


inner  product  is  defined  by 
= [ u(x)  v(x)  c~2(x)  p-1(x)dx 


(2.44) 


(Au,v)  = (u,Av) 


(2.45) 


for  all  u and  v in  the  domain  of  A which  satisfy  the  boundary 
condition  (2.9)  or  (2.10).  Moreover,  (2.44)  defines  the  inner 
product  in  the  Hilbert  space  3C  = L2 (ft,c-2 (x) p-1 (x)dx)  of  functions 
on  ft  which  are  square-integrable  with  respect  to  the  measure 
c-2(x)p_1 (x)dx.  Hence  (2.45)  implies  that  A,  acting  in  the  clas- 
sical sense  on  functions  which  satisfy  (2.9)  or  (2.10),  is  a 
symmetric,  operator  in  3f.  Moreover,  (2.42)  implies  that 


(Au,u) 


3u  3u  - 

n 3x.  3x.  P 
A 3 J 


(x)dx  > 0 


(2.46) 


for  all  u in  the  domain  of  A.  Hence  A is  positive.  It  was  shown 
in  [42]  and  [43]  how  the  domain  of  A could  be  enlarged  to  obtain 
an  extension  A of  A which  is  selfadjoint  and  positive  in  3C.  The 
boundary  condition  (2.9)  or  (2.10)  is  incorporated  into  the  defi- 
nition of  the  domain  of  A.  Moreover,  the  construction  provides  a 
meaningful  generalization  of  the  boundary  conditions  for  arbitrary 
domains  ft  C R3,  The  precise  definitions  and  results  are  reviewed 
in  sections  3-7  below. 


The  operator  (2.38)  for  acoustic  waves  in  solids  can  be 
treated  similarly.  The  divergence  theorem  implies 


ft 


(Au).  v,  p(x) dx 


'j  3 

I, 


Jim,  . 3uJt  , 

C (x)  -r- — -T-— dX  - 
0 jk  9xm  3xk 


3(1 


f-m,  . 3u 

cjk(x)  3x 


(2.47) 

dUJl  ) 

3-\l  VS 

m j J 


It  follows  that  if  an  inner  product  is  defined  by 


(u,v) 


u (x)  v (x) 

h 3 3 


p(x)dx 


(2.48) 


then  (2.45)  holds  for  all  u and  v in  the  domain  of  A which  satisfy 
the  boundary  condition  (2.20)  or  (2.21).  Moreover,  (2.48)  defines 
the  inner  product  in  the  Hilbert  space  = L2 (H ,C3 ,p(x)dx)  of 
functions  from  il  to  C3  which  are  square  integrable  with  respect  to 
the  measure  p(x)dx.  Hence  A,  acting  in  the  classical  sense  on 
functions  which  satisfy  (2.20)  or  (2.21),  is  a symmetric  operator 
in  J(.  Moreover,  (2.47)  implies  that 


I 


j 
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(Au,u)  = 


•to 


0 3u0  9u. 

x.m,  . 1 i , ^ „ 

cjk(x)  3x~  jrdx>-° 
m k 


(2.49) 


for  all  u in  the  domain  of  A by  the  assumed  positivity  of  the 
energy  density,  (2.35).  It  will  be  shown  in  section  8 below  how 
the  domain  of  A can  be  enlarged  to  obtain  a selfadjoint  positive 
extension  A of  A. 


A Hilbert  space  and  selfadjoint  positive  operator  A on  If 
can  be  associated  with  each  acoustic  wave  propagation  problem  by 
the  method  indicated  above.  A theory  of  solutions  of  the  initial- 
boundary  value  problem  (2.39)  - (2.41)  may  then  be  based  on  A in 
the  following  way.  First  of  all,  the  problem  can  be  formulated  as 
an  initial  value  problem  in  Jf.  A function  u:  P.  -*  3f  is  sought 

such  that 


+ Au  = 0 for  all  t £ R 


(2.50) 


u(0) 


= f and 


du(0) 

dt 


g in  3f 


(2.51) 


The  spectral  theorem  for  A: 

A = £ A dJI(A)  (2.52) 

and  the  associated  operator  calculus  make  it  possible  to  construct 
the  generalized  solution 

u(t)  = (cos  t A1/>2)f  + (A~1/f2  sin  t A!^2)g  (2.53) 

The  coefficient  operators  in  (2.53)  are  bounded  and  hence  u(t)  is 
defined  for  all  f and  g in  3f  and  defines  a curve  in  C(R,K),  the 
class  of  continuous  Jf-valued  functions  on  R.  The  differentiability 
properties  of  u(t)  depend  on  those  of  f and  g.  Two  cases  will  be 
mentioned. 


2.4  Solutions  in  3f 

If  f € JC  and  g € 5C  then  u(t)  is  continuous  in  7f  and  u(0)  = f. 
However,  u(t)  will  not  in  general  be  differentiable,  and  hence 
(2.50)  and  the  second  initial  condition  need  not  hold.  In  this 
case  u(t)  coincides  with  the  "generalized  solution  in  3 C which  was 
defined  and  studied  by  M.  Vishik  and  0.  A.  Ladyzhenskaya  [32]. 
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2.5  Solutions  with  finite  energy 

If  f £ D(a'/2)  and  g E Jf  then  u is  in  the  class 

C1  (R,?0  n C(R,D(A1/2))  (2.54) 

This  follows  easily  from  (2.53)  and  the  spectral  theorem.  Hence, 
u satisfies  (2.51)  but  (2.50)  need  not  hold.  In  this  case  uft) 
coincides  with  the  "solution  with  finite  energy"  which,  for  arbi- 
trary domaii.s  ft,  was  defined  and  studied  by  the  author  in  [33,34, 
42].  The  existence  and  uniqueness  of  solutions  with  finite  energy 
was  proved  in  [33,34]. 


3.  PROPAGATION  IN  HOMOGENEOUS  FLUIDS 


Propagation  in  an  unlimited  homogeneous  fluid  is  analyzed  in 
this  section.  In  the  notation  of  section  2 this  is  the  special 
case  where  ft  = R3  and  p(x)  = p and  c(x)  = c are  constant  for  all 
x 6 R3.  It  will  be  enough  to  treat  the  case  c = 1 since  the  gen- 
eral case  can  be  reduced  to  this  one  by  the  change  of  variable 
ct  -*■  t.  With  these  simplifications  the  wave  equation  (2.8)  reduces 
to  the  d'Alembert  equation 


Vu  , ?iu  9_V 

H 3xl 


= 0 


(3.1) 


and  tee  propagation  problem  is  simply  the  Cauchy  problem  for  (3.1). 
The  spectral  and  asymptotic  analysis  of  solutions  in  L2(R3)  of 
(3.1)  was  developed  in  detail  in  [42].  Only  the  principal  concepts 
and  results  are  reviewed  here. 


The  operator  in  L2(R3)  defined  by  Au  = ~(32u/9x2 + 92 u/9x2 + 
92u/9x2)  acting  in  the  domain  D(A)  = P(R3),  the  L.  Schwartz  space 
of  testing  functions,  is  known  to  be  essentially  selfadjoint  [18] 
Thus  A has  a unique  selfadjoint  extension  in  L2(R3)  which  will  be 
denoted  here  by  A0 . This  operator  may  be  defined  by 


D(A0)  = L2(R3)  n < 


u: 


92 

9x? 


u 9 2 u 

9x2 


92u 

9x| 


L2(R3) 


(3.2) 


and 


A0u  = 


9^u  32u  32u  ' 

9x2  3^5“  3^' 


for  all  u e D(A0) 


(3.3) 


where  the  derivatives  are  interpreted  in  the  sense  of  Schwartz's 
theory  of  distributions.  A0  is  known  to  be  non-negative  and  it 
is  obviously  real;  that  is 


1 
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AQu  = Aqu  for  all  u € D(Afl) 
where  the  bar  denotes  the  complex  conjugate. 


(3.4) 


The  d'Alembert  equation  (3.1)  will  be  interpreted  as  the 
equation 


d2u 

dt2 


+ Aqu  = 0 


(3.5) 


for  an  L2 (R3)-valued  function.  Hence  the  solution  in  L2(R3)  of 
the  Cauchy  problem  can  be  written 

u(t)  = (cos  t Ag/2)f  + (A“1/2sin  t A0V2)g  (3.6) 

where  u(0)  = f and  du(0)/dt  = g are  in  L2(R3).  If  it  is  assumed 
that  f(x)  and  g(x)  are  real-valued  and 

f e L2(R3),  g G D(A^1/2)  (3.7) 

then  it  follows  from  (3.4)  that 

u(t,x)  = Re  (v(t,x)}  (3.8) 

where 

v(t,*)  = exp  (-itA^2)h  (3.9) 

and 

h - f + i A'1'2  g G L2(R3)  (3.10) 

In  what  follows  attention  is  restricted  to  this  case. 


An  eigenfunction  expansion  for  A0  may  be  based  on  the 
Plancherel  theory  of  the  Fourier  transform  in  L2(R3).  If 


wQ(x,p)  = 


(2t0  3/2 


exp  (i  x • p) , p G R3 


(3.11) 


where  x • p = x3Pj  + x2p2  + x3p3  then  the  main  results  of  the  theory 
state  that  for  all  f G L2(R3)  the  following  limits  exist 


f(p)  2 (*0f)(p)  = L2(R3)-lim 


M-KO 


f (x)  = ($;f)(x)  = L2(R3)-lim 


w0(x,p)  f (x)dx 


x <M 


> (3.12) 


w.(x,p)  f (p) dp 


IPl<M 
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and  <J>0:  L2(k3)  -+•  L2(R3)  is  unitary.  These  relations  will  often 

be  written  in  the  symbolic  form 


■ f(p) 


w0(x,p)  f(x)dx,  f(x) 

JR3 


, w0(x,p) 
JR3 


f(p)dp  (3.13) 


but  must  be  interpreted  in  the  sense  (3.12).  The  utility  of  the 
Fourier  transform  is  due  to  the  fact  that  if  f and  9f/9x.  are  in 
L2(R3)  then  ^ 

*o  = iPj  f(P)»  3 = 1.2,3  (3.14) 

In  particular,  it  follows  that 

D(A0)  = L2(R3)  n {u:  |p|2u(p)  € L2(R3)}  (3.15) 

Aq  has  the  spectral  representation 

r 

A0  = Jp  X dllp(A)  (3.16) 

with  spectral  family  {no(A)}  defined  by 

no(A)f(r.)  = f Wp  (x, p)  f (p) dp,  A > 0 (3.17) 

'|p|<A 

It  follows  that  A0  is  an  absolutely  continuous  operator  [18,42] 
whose  spectrum  is  the  interval  [0,ro). 

The  above  results  imply  that  <J>0  defines  a spectral  represen- 
tation for  Afl  and  functions  of  AQ . In  particular,  if  T(A)  is  any 
bounded  Lebesgue-measurable  function  of  A > 0 then 

$04'(Ap)f(p)  = T(  | p j ?)  f (p)  (3.18) 

These  results  imply  that  the  wave  function  v(t,x)  defined  by 
(3.9)  has  the  representation 

v(t,x)  = [ Wp (x , p)  exp  (-it  j p | ) h (p) dp  (3.19) 

JR3 

The  function  wQ(x,p)  is  a generalized  eigenfunction  for  A0 . This 
means  that  w0(*,p)  is  locally  in  D(Aq);  i.e.,  <J>w0(*,p)  G D(A0)  for 
every  4>  6 P(R3)  and 

Aowo<*»P>  = A wo^*.P^>  A = I p 1 2 (3.20) 


The  functions 
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w0(x,p)  exp  (-it  | p | ) = -(-2-^2  exp  {i(x*  p - t|p|)}  (3.21) 

are  solutions  of  the  d'Alembert  equation  which  represent  plane 
waves  propagating  in  the  direction  of  the  vector  p £ R3.  Hence, 
(3.19)  is  a representation  of  a localized  acoustic  wave  as  a 
superposition  of  the  elementary  waves  (3.21). 

The  spectral  integral  (3.19)  is  the  starting  point  for  the 
asymptotic  analysis  of  the  behavior  for  t ■*  00  of  solutions  in 
L2(R3)  of  the  d'Alembert  equation.  It  is  convenient  to  begin  the 
analysis  with  the  special  case  where  h is  in  the  class 

VQ  (R3)  = P(R3)  n (h:  h(p)  = 0 for  |p|  < a,  a=  a(h)  > 0}  (3.22) 

The  analysis  will  then  be  extended  to  the  general  case  by  using 
the  easily  verified  fact  that  P0(R3)  is  dense  in  L2(R3). 

A , A 

If  h 6 Vq(R3)  and  the  support  of  h satisfies 

6upp  h C {p:  0 < a < |pj  < b}  (3.23) 

then  the  spectral  integral  (3.19)  converges  both  in  L2(R3)  and 
pointwise  to  v(t,x)  and 


v(t,x) 


1 

(27P77 


exp  { i (x  • p - t[p|))h(p)dp  (3.24) 


' a< | p | <b 

To  find  the  behavior  of  v(t,*)  £ L2(R3)  for  t 00  introduce 
spherical  coordinates  for  p: 


p = pw,  p > 0,  w £ S2,  dp  * p2dpdco 


(3.25) 


where  S2  represents  the  unit  sphere  in  R3  with  center  at  the  origin 
and  dtu  is  the  element  of  area  on  S2.  This  gives  the  representation 


v(t,x) 


1 

( 2tt)  3/2 


b 

e"ltP  V(x , p)  p2dp 
a 


(3.26) 


where 


V(x,p) 


ipx*u 

e 


A 

h(pa)du) 


V 


(3.27) 


The  asymptotic  behavior  of  V(x,p)  for  |x|  ■*  00  will  be  calculated  and 
used  to  find  the  behavior  of  v(t,x)  for  t •+  °°.  Application  of  the 
method  of  stationary  phase  [2,23]  to  (3.27)  with  x = rO,  r > 0, 

8 £ S2  implies  that  if 
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V(x,p)  = 


2u 


iprj 

+ q„(x,p) 


cipr  h(P0)  + 


2rr 

-iprj 


e_iPr  h(-P0) 


(3.28) 


then  there  exists  a constant  M0  = Mfl(h)  such  that 

J q0  Cx,p)  | < MQ/r2  foi  all  r > 0,  a < p < b and  0 g S2  (3.29) 

Substituting  (3.28)  into  (3.26)  gives 

v(t,x)  = G(r  - t,0)/r  + G’(r+  t,0)/r  + q^t.x)  (3.30) 

where  G(t,0)  and  G'(t,0)  are  the  functions  of  T 6 R and  0 6 S2 
defined  by 


G(t,0) 


and 


1 

(2tt)  1/2 


b 

elTp  h(P0) (-ip)dp 
a 


(3.31) 


G'(t,0) 


1 

( 2tt)  1/2 


■ —a 

elTp  h(p0) (-ip) dp 


(3.32) 


J -b 

Moreover,  the  estimate  (3.29)  implies  that  qj/t.x)  satisfies 

I qj ( t , x) | < Mj/r2  for  all  r > 0,  t E R and  0 e S2  (3.33) 
where  Mt  = M,  (h)  = (2tt)~3/2  (b3  - a3)  M0(h)/3. 

The  principal  result  of  this  section  states  that 
v (t,x)  = G(r-  t,0)/r,  x = r0  (3.34) 

is  an  asymptotic  wave  function  for  v(t,x)  in  L2(R3);  that  is, 

lim  II  v(t,  •)  - v°°( t , 4 ) II T = 0 (3.35) 

t-HX3  L2^K  / 


Before  indicating  a proof  it  is  necessary  to  complete  the  statement 
of  the  theorem  by  defining  the  profile  G for  arbitrary  h e.  L„(R3). 
When  h £ P0(R3),  G is  defined  by  (3.31)  and  a simple  calculation 
gives 


I Gil 


2 

L. 


(RxS2) 


J a< 


_ I P 


h (p)  | 2dp  = I hi 
<b 


l2(r3) 


= I hi 


L2(R3) 


(3.36) 


r 


r 
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Hence  the  correspondence 

h -*■  G = 0h  G L2(R  x S2)  (3.37) 

can  be  extended  to  all  h e L2(R3)  by  completion.  Another  method, 
based  on  the  Plancherel  theory  in  L2(R,L2(S2))  is  given  in  [42]. 
It  is  not  difficult  to  verify  by  constructing  0-1  that 

0:  L2(R3)  -*•  L2(R  x S2)  is  unitary  (3.38) 

A similar  extension  of  the  definition  (3.32)  of  G'  may  be  made. 

A proof  of  (3.35)  will  new  be  outlined.  Note  first  that  the 
function  G'(r+t,0)/r  tends  to  zero  in  L2(R3)  when  t 00 . This 

follows  from  the  simple  calculation 


| G * (r+  t,0)/r  j 2dx  = 
R3 


|G'(r+  t,0)  | 2d0dr 

o is2 


|G'(r,0) |2d0dr 

J t 's2 


(3.39) 


and  the  fact  that  G'  G L2(R  x S2) . The  proof  that,  in  (3.30), 
qj  (t,  •)  + 0 in  L2(R3)  when  t ->■  » is  based  on  the  following  lemma. 


3.1  Convergence  lemma 

Let  ft  C R3  be  an  unbounded  domain  and  let  u(t,x)  have  the 
properties 

u( t , * ) G L2 (ft)  for  every  t > t0  (3.40) 

lim  lu(t,OllL  /K,  = 0 for  every  compact  K C R3  (3.41) 

t -K»  ' 

|u(t,x)|  < M/ | x | 2 for  every  |x|  > rQ  (3.42) 

where  t0 , r0  and  M are  constants.  Then 

lim  Hu(t,*)ll^  = 0 (3.43) 

t-x»  ' 


Only  the  case  ft  = R3  of  the  lemma  is  needed  here.  The  more 
general  case  is  used  in  section  4.  A simple  proof  of  the  lemma 
is  given  in  [42]. 

A _ 

The  proof  of  (3.35)  for  the  case  h G #0(R3)  may  be  completed 
by  applying  the  lemma  to  u(t,x)  = q;(t,x).  (3.33)  states  that  qj 


1 


i 
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satisfies  (3.42)  while  (3.40)  and  (3.41)  follow  from  (3.30).  To 
verify  (3.41)  note  that  G'(r+t,6)/r  satisfies  it  by  (3.39). 
Moreover,  if  K C {x:  |x|  < R)  then  by  direct  calculation 


| |G(t 

* 

i 

,1 


- t,0)/r | 2dx  < 


,0)/r  2dx 


|G(r-t, 

J |x | <R 

r R-u 

t,0) | 2d0dr  = |G(r,0) |2d0dr 

J -t  Js2 


(3.44) 


The  last  integral  tends  to  zero  when  t because  G € L2(R  * S2)  . 
Finally,  v(t,x)  satisfies  (3.41).  When  h £ D0(R3)  this  can  be 
verified  directly  from  (3.24)  by  an  integration  by  parts. 

A 

The  proof  of  (3.35)  indicated  above  is  valid  when  h £ DQ(R3). 

To  prove  (3.35)  for  general  h £ L2(R3)  note  that 

v(t,*)  = U0(t)h  where  UQ(t)  = exp  (~itA01/2)  (3.45) 

is  unitary.  In  particular, 

B UQ (t) B = 1 for  all  t £ R (3.46) 

Similarly,  if  UQ(t):  L2(R3)  -*  L2(R3)  is  defined  by 

v°°  ( t , • ) = U^(t)h  (3.47) 

CO 

then  it  follows  from  (3.44)  and  (3.36)  that  UQ(t)  is  contractive: 

iu7(t)l  < 1 for  all  t £ R (3.48) 

The  general  case  of  (3.35)  now  follows  from  the  special  case 
h £ V (R3),  the  density  of  D0(R3)  in  L2(R3)  and  the  estimates 
(3.46)  and  (3.48).  The  details  are  given  in  [42]. 

The  real  part  of  the  asymptotic  wave  function  (3.34)  is  another 
function  of  the  same  form.  Hence,  (3.8)  and  (3.35)  imply  a similar 
result  for  the  solution  in  L2(R3)  of  the  Cauchy  problem.  The  result 
may  be  formulated  as  follows. 


3.2  Theorem 

Let  f and  g be  real -valued  functions  such  that  f £ L2(R3)  and 
g £ D(Ag  1/,z)  . Let  u(t,x)  be  the  corresponding  solution  in  L2(R3) 
of  the  d'Alembert  equation  given  by  (3.6).  Define  the  asymptotic 
wave  function 


I 

r 

j 

\i 
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u (t,x)  = — 


F(r-  t,8) 


x = r0 


(3.49) 


where 


F(t,0)  = Re  {G(t,0) } 


G = 0h  = 0(f  + iA71/2g) 


(3.5J ) 


lim  II u ( t , • ) - u (t,*)ll,  /r. 3 \ 

t-x»  “2'  ' 


(3.52) 


Stronger  forms  of  convergence  than  (3.52)  can  also  be  proved 
under  suitable  hypotheses  on  the  initial  state.  In  particular, 
convergence  in  energy  holds  if  the  initial  state  has  finite  energy. 
A result  of  this  type  is  formulated  at  the  end  of  section  4 for 
the  more  general  case  of  an  initial-boundary  value  problem  for  the 
d'Alembert  equation  in  an  exterior  domain. 


4.  SCATTERING  BY  OBSTACLES  IN  HOMOGENEOUS  FLUIDS 


The  scattering  of  localized  acoustic  waves  by  bounded  rigid 
obstacles  immersed  in  an  unlimited  homogeneous  fluid  is  analyzed 
in  this  section.  The  corresponding  boundary  value  problem  is 


92u  f 92u  a2u  32u 

9t2  i 9x^  9x|  9x| 


0 for  t > 0,  x € JJ 


(4.1) 


= 0 for  t > 0,  x £ 9(2 


(4.2) 


u(0,x)  = f(x)  and  9u(0,x)/3t  = g(x)  for  x £ (2 


(4.3) 


where  (2  C R3  is  an  exterior  domain  (i.e.,  T = R3  - (2  is  compact). 
This  problem  will  be  treated  as  a perturbation  of  the  Cauchy  prob- 
lem of  section  3. 


A formulation  of  the  initial-boundary  value  problem  (4.1)  - 
(4.3)  which  is  applicable  to  arbitrary  domains  (2  C R3  was  given  by 
the  author  in  [33,42].  That  work  provides  the  starting  point  for 
the  analysis  of  this  section  and  sections  5 and  6.  The  principal 
definitions  and  results  are  summarized  here  briefly. 

The  formulation  makes  use  of  the  Hilbert  space  L2((2)  and  the 
following  subsets  of  L2(lf2). 


T 
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L*(ft)  = L2(fi)  n {u,  3u/9Xj  e L2(fl)  for  j = 1,2,3}  (4.4) 

L2(A,(2)  = L2(fi)  n (u:  Au  e L2((2)}  (4.5) 

L£(A,(l)  = L*(S2)  n L2(A,«)  (4.6) 

where  Au  = 92u/9x2  + 92u/9x2  + 92u/3x2  denotes  the  Laplacian  of  u. 
The  derivatives  in  these  definitions  are  to  be  interpreted  in  the 
sense  of  the  theory  of  distributions.  The  sets  (4.4),  (4.5)  and 
(4.6)  are  linear  subsets  of  L2(ft).  Moreover,  they  are  Hilbert 
spaces  with  inner  products 

3 

(u,v)j  = (u,v)  + l (3u/9Xj ,9v/3xj)  (4.7) 

j=l 

(u,v)A  = (u , v)  + (Au, Av)  (4.8) 

(u , v)  . = (u , v)  j + (Au,Av)  (4.9) 

respectively,  where  (u,v)  is  the  inner  product  in  L2(i2). 


4.1  Definition 

A function  uS  L2(A,ft)  is  said  to  satisfy  the  generalized 
Neumann  condition  for  Q if  and  only  if 

3 

(Au ,v)  + 7 (3u/9x . , 9v/3x . ) = 0 for  all  v£  L.i((2)  (4.10) 

L J 1 

J=1 

Note  that  (4.10)  defines  a closed  subspace 

L^(A,f2)  = L2(A,fi)  n (u:  u satisfies  (4.10)}  (4.11) 

i N 

in  the  Hilbert  space  L*(A,(2).  The  condition  "u  € L2(A,Q)"  is  a 
generalization  of  the  Neumann  boundary  condition  (4.2).  It  is 
meaningful  for  arbitrary  domains  12.  Moreover,  it  reduces  to  (4.2) 
whenever  is  sufficiently  smooth  (see  [42,p.41]  for  a discussion). 

The  construction  of  solutions  of  the  initial-boundary  value 
problem  (4.1)  - (4.3)  given  below  is  based  on  the  linear  operator 
A = A(I2)  in  L2(ft)  defined  by 

D(A)  = lJ(A,«)  (4.12) 

Au  = -Au  for  all  u e D(A)  (4.13) 

The  utility  of  this  operator  is  based  on  the  following  theorem 
which  is  proved  in  [42]. 
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4.2  Theorem 


A Is  a selfadjoint  real  positive  operator  in  L2(ft). 
D(Al/z)  = L*(ft)  and 

3 

0 A 1/2 ull  2 = l II  8u/8x J 2 for  ail  u€  D(A1/2) 

j=l 


Moreover , 


(4.14) 


The  operator  A may  be  used  to  construct  "solutions  in  L2(ft)" 
and  "solutions  with  finite  energy"  of  (4.1)  - (4.3),  as  described 
in  section  2.  The  solution  in  L2  (ft)  will  be  considered  here.  As 
in  section  3,  if  f 6 L2 (ft)  and  g G D(A-1/2)  then 


u(t,x)  = Re  { v(t ,x) } 


(4.15) 


where 

v ( t , • ) = exp  (-itA1/2)h,  h = f + iA~1/2g  (4.16) 

The  properties  of  the  operator  A stated  in  the  theorem  above 
are  valid  for  arbitrary  domains  ft  C R3.  It  was  shown  in  [42]  that 
if  ft  is  an  exterior  domain  then  A has  a continuous  spectrum. 
Moreover,  if  ft  has  the  local  compactness  property  (defined  below) 
then  there  exist  eigenfunction  expansions  for  A in  terms  of  gener- 
alized eigenfunctions  which  are  perturbations  of  the  plane  wave 
eigenfunctions  of  section  3.  In  the  remainder  of  this  section  the 
eigenfunction  expansions  are  described  and  used  to  analyze  the 
structure  of  solutions  of  the  scattering  problem  (4.1)  - (4.3). 

The  principal  result  of  the  analysis  states  that  the  behavior  of 
the  acoustic  field  for  large  times  is  described  by  an  asymptotic 
wave  function  of  exactly  the  same  form  (3.49)  as  when  there  is  no 
obstacle.  The  only  effect  of  an  obstacle  is  to  modify  the  wave 
profile  F(x,0).  Moreover,  a procedure  is  given  for  calculating 
the  modified  profile  when  the  obstacle  and  the  initial  state  are 
known. 


4.3  Distorted  plane  wave  eigenfunctions 

Two  families  of  generalized  eigenfunctions  of  A,  denoted  by 
w+(x,p)  and  w_(x,p)  respectively,  were  defined  in  [42].  They  are 
perturbations  of  the  plane  wave  eigenfunctions  w„(x,p)  and  have 
the  form 

w±(x,p)  = wQ(x,p).  + w|(x,p),  p G R3  (4.17) 

where  w.J_(x,p)  and  w^(x,p)  may  be  interpreted  as  secondary  fields 
which  are  produced  when  the  obstacle  T = R3  - ft  is  irradiated  by 
the  plane  wave  wQ(x,p).  Mathematically,  w+(x,p)  and  w_(x,p)  must 
satisfy 
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(A  + | p j 2)  w±(x,p)  = 0 for  x e J2  (4.18) 

9w+(x,p) 

= 0 for  xe  3!)  (4.19) 

9v 

However,  they  are  not  completely  determined  by  these  conditions. 
Instead,  w+(x,p)  is  determined  by  (4.18),  (4.19)  and  the  condition 
that  w_|.(x,p)  should  describe  an  outgoing  secondary  wave.  This  is 
implied  by  the  Sommerfeld  condition  for  outgoing  waves: 

9w!(x,p) 

i I P I w|(x,p)  = 0 ( | x | - 1 ) , | x | -*•  oo 

8|x| 

w|(x,p)  = 0 ( | x | ~ 1 ) , j x | -*■  “ 

Similarly,  w_(x,p)  is  determined  by  (4.13),  (4.19)  and  the  condi- 
tion that  w^(x,p)  should  describe  an  incoming  secondary  wave, 
which  is  implied  by  the  Sommerfeld  condition  for  incoming  waves: 

9w^(x,p) 

+ i | p | w^(x,p)  = 0 ( | x | *),  | x | -*  o° 

w^(x,p)  = 0(  | x | _ 1 ) , | x | ->  00 

Of  course,  if  9 ft  is  not  smooth  then  the  boundary  condition 

(4.19)  must  be  understood  in  the  generalized  sense  of  (4.10).  A 
technical  difficulty  is  caused  by  the  fact  that  w+(*,p)  cannot  be 
in  D(A)  = Lj(A,ii)  because  the  spectrum  of  A is  continuous.  This 
is  overcome  by  requiring  that 

<K(*,p)  6 I^(A.fl)  (4.22) 

for  all  <{)  6 D(R3)  such  that  4>(x)  = 1 in  a neighborhood  of  9fl. 
Generalized  eigenfunctions  with  these  properties  will  be  called 
"distorted  plane  waves,"  following  T.  Ikebe  [16]. 

The  uniqueness  of  distorted  plane  waves  satisfying  (4.18), 

(4.20)  or  (4.21)  and  (4.22)  was  proved  in  [42]  for  arbitrary  ex- 
terior domains.  However,  to  prove  their  existence  it  was  necessary 
to  impose  a condition  on  9f2.  To  define  it  let 

11R  ■ 11  n Cx:  | x | < R]  (4.23) 

L^°c(f2)  = {u:  u € L203r)  for  every  R > 0} 


(4.24) 


24 


L*,£oC(ft)  = L£°C(ft)  n {u:  au/axj  G L^°C(ft)  (4.25) 

for  j = 1,2,3) 


and  define  the 


4.4  Local  compactness  property 

A domain  ft  C R3  is  said  to  have  the  local  compactness  property 
if  and  only  if  for  each  set  S C L1’^oc(ft)  and  each  R > 0 the 
condition 


< C(R)  for  all  u G S 


(4.26) 


implies  that  S is  precompact  in  L2  (ft^) ; i.e.,  every  sequence  {u^} 
in  S which  satisfies  (4.26)  has  a subsequence  which  converges  in 
L2(ftft).  The  class  of  domains  with  the  local  compactness  property 
will  be  denoted  by  LC. 

The  local  compactness  property  is  known  to  hold  for  large 
classes  of  domains.  S.  Agmon  has  proved  it  for  domains  with  the 
"segment  property"  [1] . A generalization  of  the  segment  property, 
called  the  "finite  tiling  property"  was  given  by  the  author  in 
[42].  As  an  application  of  this  condition  it  can  be  shown  that 
the  local  compactness  property  holds  for  the  many  simple,  but  non- 
smooth, boundaries  that  arise  in  applications,  such  as  polyhodra, 
finite  sections  of  cylinders,  cones,  spheres,  disks,  etc.  The 
following  existence  theorem  was  proved  in  [42] . 


4.5  Theorem 

Let  ft  C R3  be  an  exterior  domain  such  that  ft  G LC . Then  for 
each  p G R3  there  exists  a unique  outgoing  distorted  plane  wave 
w+(x,p)  and  a unique  incoming  distorted  plane  wave  w_(x,p) . 

The  outgoing  (resp.  incoming)  property  of  v’(x,p)  (resp. 
w'(x,p))  is  made  explicit  by  the  following  corollary. 


4.6  Corollary 

Under  the  same  hypotheses  there  exist  functions 
T+(Q,p)  6 C (S2  x {R3  - 0})  such  that 


w|(x,p)  = 


r 


T+(0,p)  + w"(x,p) , x = r0 


(4.27) 
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whe  re 


w"(x,p)  = 0(r-2) , r -*•  » 


(4.28) 


uniformly  for  0 = x/r  G S2  and  p in  any  compact  subset  of  R3  - {0} . 


In  acoustics  the  functions  T+(0,p)  and  T_(0,p)  are  called  the 
far-field  amplitudes  of  the  distorted  plane  waves. 


4.7  The  eigenfunction  expansion  theorem 


Each  of  the  families  (w+(*,p):  p G R3}  and  (w_(*,p):  p G R3} 

defines  a complete  set  of  generalized  eigenfunctions  of  A in  the 
sense  described  by  the  following  theorems. 


4.8  Theorem 

For  each  f G L2(Q)  the  following  limits  exist 


f±(p)  = L2(R3)-lim 

M-*°° 


a 


w+(x,p)  f (x) dx 


M 


(4.29) 


f(x)  = L.(ft)-lim 


M-*°  J | p | <m  ± 


w+(x,p)  f±(p)dp 


where  fyj  = n (x:  |x|  < M) . Moreover,  the  operators 

$>+:  L?((2)  L2(R3)  defined  by 


$+f  = f+ 


(4.30) 


are  unitary. 

The  relations  (4.29)  will  usually  be  written  in  the  symbolic 


form 


r 


w+(x, p)  f (x)dx,  f(x)  = , w+(x,p)  f+ (p) dp 

0 ‘ •'R3  " ” (4.31) 


f+(p)  = 

but  must  be  understood  in  the  sense  of  (4.29). 


4.9  Theorem 

If  {11(A)}  denotes  the  spectral  family  of  A: 


A dJI(A) 


(4.32) 


J 
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then  TI (X)  has  the  eigerf unction  expansions 


n<A)  f(x> 


W±(x,p)  f+ (p)dp,  X > 0 

■ |p|<^ 


(4.33) 


In  particular,  A is  an  absolutely  continuous  operator  whose  spectrum 
is  the  interval  (0,°°). 

The  last  result  implies  that  <i>+  and  4>_  define  spectral  repre- 
sentations for  A in  the  sense  of  the  following  corollary. 


4.10  Corollary 

If  'i'(A)  is  a bounded  Lebesgue-measurable  function  of  X > 0 
then  for  all  f 6 L2 (ft) 

*±¥(A)f(p)  = V(|p|2)  f±(p)  (4.34) 

These  results  provide  a complete  generalization  of  the 
Plancherel  theory  to  exterior  domains  ft  G LC. 


4.11  The  eigenfunction  expansions  and  scattering  theory 

The  results  stated  above  imply  that  the  wave  functions 

v(t,*)  = exp  (-itA1/,2)h,  he  L2(ft)  (4.35) 

have  the  spectral  integral  representations 

v(t,x)  = [ w+(x,p)  exp  (-it  j p I ) h+(p)dp  (4.36) 

1R3  " 

Note  that  (4.36)  defines  two  representations,  corresponding  to 
w+(x,p)  and  w_(x,p).  They  will  be  called  the  outgoing  and  incoming 
representations,  respectively. 

The  representations  (4.36)  and  the  results  of  section  3 will 
now  be  used  to  derive  the  asymptotic  behavior  of  v(t,x)  for  t -*  00 . 

To  begin  consider  an  initial  state  h € L2(ft)  such  that 

h_  G VQ (R3)  (4.37) 

Such  states  are  dense  in  L2 (ft)  because  0g(R3)  is  dense  in  L2(R3) 
and  L2(ft)  -*■  L2(R3)  is  unitary.  The  wave  function  corresponding 

to  (4.37)  is 
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r(t,x)  = w (x,p)  exp  (-it|p|)  h 

JR3 


where  the  integral  converges  both  pointwise  and  in  L2(fi)  to  v(t,x). 
To  discover  the  behavior  of  v(t,x)  for  t -*■  °°  substitute  the 
decompositions  (4.17)  and  (4.27)  for  w_(x,p)  into  (4.38)  and  write 


v(t,x)  = v0(t,x)  + v'(t,x)  + v" ( t , x) 


where 


:,x)  = 

R 


w (x,p)  exp  (-it | p I ) h (p)dp 


(4.39) 


(4.40) 


■y 

v’(t,x)  =-  exp  t-i | p | (r  + t) } T (0,p)  h (p)dp  (4.41) 

•'R3 


/’(t,x)  = w" (x > p)  exp  (—it | p ( ) h 

R3 


(p)  dp 


(4.42) 


Note  that  vfl(t,x)  is  a solution  in  L2(R3)  of  the  d'Alembert  equa- 
tion. Indeed,  h_  = $_h  = $0  (^^.h)  = hfl  where 


hfl  = 4>*$>_h  e L2(R3) 


(4.43) 


v0(t,*)  = w0(*,p)  exp  (-it | p | ) h. (p) dp 

JR3 


(4.44) 


= exp  (-itA3/2)h0 


Thus  vQ(t,x)  represents  a wave  in  an  unlimited  fluid  containing  no 
obstacles.  It  will  be  shown  that  v(t,x)  is  asymptotically  equal 
to  this  wave  when  t •+  m;  i.e., 


lim  « v(t, -)  - v0(t,*)lL2(fl)  = 0 


(4.45) 


To  see  this  note  that,  in  (4.39),  v'(t,x)  has  the  form 
v ' ( t , x)  = G' (r  + t , 0 ) / r 


(4.46) 


It  was  shown  in  section  3 that  such  functions  tend  to  zero  in 
L2(R3)  when  t -*■  00  (see  (3.39)).  It  is  easy  to  check  that  (4.37) 
implies  that  G'  G L2(R  * S2) . Finally,  condition  (4.28)  for 
w"(x,p)  implies  that  the  term  v"(t,x)  in  (4.39)  satisfies 
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|v"(t,x) | < M/ | x | 2 for  all  |x|  > 0 and  t £ R (4.47) 

with  a suitable  constant  M.  Hence,  the  convergence  lemma  of 
section  3,  applied  to  v"  = v - v0  - v'  implies  (4.45)  if  v"(t,x) 
satisfies  the  local  decay  condition  (3.41).  For  v'(t,x)  this 
condition  follows  from  (4.46).  For  v(t,x)  and  vQ(t,x)  it  follows 
from  the  local  compactness  property.  A proof  may  be  found  in 
[42].  Thus  (4.45)  is  established  for  all  h_  £ D0(R3).  The  main 
result  of  this  section  is  the 


4.12  Theorem 

For  all  h e L2(&)  if  v(t,*)  = exp  (-itA^2)h  and  vo(t,») 
= exp  (-itAg/2)  ($*4>_)h  then 


lim  |v(t,-)  - v0(t,.)lL2(ft)  = 0 


(4.48) 


This  result  follows  immediately  from  the  special  case  (4.37) 
proved  above,  the  density  of  P0(R3)  in  L2(R;>)  and  the  unitarity 
of  the  operators  exp  (-ltA1/2),  exp  (-itAg'2),  <J>0  and 

4.13  Corollary 

If  L2(ft)  L2(R3)  is  defined  by  Jqu(x)  = u(x)  for  all 

x £ £)  and  J^u(x)  - 0 for  all  x £ R3  - ft  then  the  strong  limit 

W+  = W+(A3/2 ,A3/2, J^)  = s-lim  exp  (itA^2)  exp  (-itA1^2)  (4.49) 


exists  in  L2  (fi)  and  W : L2(fi)  L2(R3)  is  given  by 


w+  = 


(4.50) 


In  particular,  W+  is  unitary. 

The  operator  W+  is  the  wave  operator  for  the  pair  aJ^2,  A1^2 
in  the  sense  of  the  time  dependent  theory  of  scattering.  The 
equivalence  of  (4.48)  and  (4.50)  is  proved  in  [42]. 


4.14  Asymptotic  wave  functions  in  L2 (ft) 


The  wave  function  in  L2(R3)  defined  by 
v q ( t , * ) = exp  (-itAg/2)h0,  h = <t>*<!>_h 


(4.51) 


has  an  asymptotic  wave  function  in  L2(R3),  by  the  results  of 
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section  3;  i.e., 

lim  II  v0  (t , * ) - v°°(t,  *)l!  , = 0 (4.52) 

t-K»  ' 


where 


v (t,x)  = G(r-t,9)/r,  x = r0 


(4.53) 


and 


G = 0hQ  = 0$*<t_h 


(4.54) 


Equations  (4.48),  (4.52)  and  the  triangle  inequality  imply  the 


4.15  Theorem 

For  each  h 6 L2(fi)  the  wave  function  v (t,*)  defined  by 
(4.53),  (4.54)  is  an  asymptotic  wave  function  in  L2  (Q)  for  v(t,*) 
= exp  (-itA1'2)h;  that  is, 

lim  U v ( t , • ) - v°°(t,*)HT  = 0 (4.55) 

t-wo 


4.16  Corollary 


The  profile  of  the 

«T,e)  - 


asymptotic  wave  function  is  given  by 

r 

e1Tp  h_(P0) (-ip)dp  (4.56) 


Jo 

where  the  integral  converges  in  I.2(R  x S2). 


This  follows  immediately  from  (4.54)  and  (3.31).  Note  that 
the  only  difference  between  the  asymptotic  wave^ functions  for  R3 
and  those  for  ft  is  that  h = <5>0h  is  replaced  by  h_  =■  h. 


4.17  Asymptotic  energy  distributions 

If  the  initial  state  h G L2(ft)  has  derivatives  in  L;(ft)  then 
the  corresponding  profile  G and  asymptotic  wave  function  v (t,x) 
will  have  corresponding  derivatives.  In  particular,  the  following 
result  was  proved  in  [42], 
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4.18  Corollary 


If  9h(x)/9x.  € l2(J2)  for  j = 1,2,3  then  3v(t,x)/9t  and 
9v(t,x)/9xj  are  in  L2 (ft)  for  all  t € R and  j = 1,2,3  and 


11m  II  9v(t,  *)/9t  - vQ  ( t , * ) II  = 0 

t-X»  2 ' ' 


lira  fl  9v(t, *)/9x. 


Vj(t'-),L!(S!> 


= 0,  j = 1,2,3 


(4.57) 


where 


V 


vk(t,x)  = Gk(r  - t,0)/r,  k = 0,1, 2, 3 

G0 (t ,0)  = -9G(t,0)/9t 

Gj (t ,6)  = -G0 (x »0)0j , 3 = 1,2,3 


(4.58) 

(4.59) 

(4.60) 


and  G(t,6)  is  given  by  (4.56). 

The  energy  integral  for  a homogeneous  fluid  is  given  by 


E(u,K,t)  = j f < 
4- 


2 3 

z 

3u(t,x) 

+ V 

3u(t,x)  1 

St 

+ L 

9x. 

3-1 

J J 

J 

> dx  (4.61) 


if  p - 1,  c = 1.  The  last  corollary  implies  that  if  u(t,x) 

= Re  (v(::,x)}  is  a solution  with  finite  energy  in  ft  then  the 
energy  in  any  measurable  cone 


C = (x  = r0:  r > 0,  0 6 C,  C S2} 


(4.62) 


has  a limit  as  t + * which  can  be  calculated  from  the  initial 
state  u(0,x)  = f(x),  9u(0,x)/3t  = g(x).  The  following  result  was 
proved  in  [42] . 


4.19  Theorem 

If  f 6 L*(ft),  g 6 L2(ft)  and  if  C is  any  measurable  cone  in 
R3  then 

lim  E(u,C  n ft,t)  = j [ | |p|f  (p)  + ig  (p) j2dp  (4.63) 

t-x»  Jc 


5.  PROPAGATION  IN  UNIFORM  TUBULAR  WAVEGUIDES 

The  propagation  and  scattering  of  localized  acoustic  waves 
is  simple  and  compound  tubular  waveguides  with  rigid  walls,  and 
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filled  with  a homogeneous  fluid,  is  analyzed  in  this  section  and 
the  next.  The  simplest  case  is  the  uniform  semi-infinite  cylinder, 
closed  by  a plane  wall  perpendicular  to  the  axis.  Other  special 


Figure  1.  Uniform  semi- inf inite  cylindrical  waveguide. 

cases  which  are  of  interest  in  applied  acoustics  include  the  cyl- 
indrical waveguide  terminated  by  a resonator,  the  tubular 


Figure  2.  Cylindrical  waveguide  terminated  by  a resonator. 

waveguide  with  a bend,  or  elbow,  coupled  cylindrical  waveguides 
with  different  cross-sections,  the  T-joint  in  a waveguide,  uniform 
waveguides  containing  an  iris,  waveguides  containing  obstacles, 
and  many  others. 

The  most  general  compound  tubular  waveguide  considered  here 
is  described  by  a domain  Q C R3  of  the  form 

!1  = 2.  u S.  U S,  u •••  u S (5.1) 

u 1 z m 


: 


Figure  6.  Waveguide  with  iris. 


where  is  a bounded  domain  and  Sj,S?, *'‘,S  are  disjoint  uniform 
semi-infinite  cylinders.  If  il  is  a waveguide,  with  rigid  walls, 
filled  with  a homogeneous  fluid,  the  corresponding  boundary  value 
problem  is  again  problem  (4.1)  - (4.3),  but  for  a domain  with  the 
structure  (5.1).  Hence,  the  Hilbert  space  formulation  of  (4.1)  - 
(4.3)  given  at  the  beginning  of  section  4,  which  is  valid  for 
arbitrary  domains  f!  C R3,  provides  a starting  point  for  the  anal- 
ysis of  the  waveguide  problems.  The  remainder  of  this  section 
presents  the  spectral  and  asymptotic  analysis  of  acoustic  waves 
in  a uniform  semi-infinite  cylindrical  waveguide.  The  general 
case  (5.1)  is  analyzed  in  section  6. 


5.1  The  uniform  semi-infinite  cylinder 

It  will  be  convenient  to  use  coordinates 
(xj ,x2 ,y)  5 (x,y)  6 R3 


(5.2) 


such  that  the  y-axis  lies  in  the  waveguide.  With  this  choice  the 
waveguide  may  be  described  by  a domain  of  the  form 


S = {(x,y):  x € G and  y > 0} 


(5.3) 


where  G C R2  defines  the  waveguide  cross  section.  It  will  be 
assumed  that  G is  bounded  and  that  S £ LC . 


The  spectral  analysis  of  the  operator  A ■=  A(S),  acting  in 
L2(S),  will  be  based  on  the  spectral  analysis  of  A(G)  acting  in 
L2 (G) . It  can  be  shown  that  the  hypothesis  S £ LC  implies  that 
G 6 LC  as  a domain  in  R2 . This  property  and  the  boundedness  of 
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G imply  that  A(G)  has  a discrete  spectrum  with  eigenvalues 

0 = XQ  < Aj  < X2  < •••  (5.4) 

such  that 


lim  X.  = 00  (5.5) 

j-*x>  J 


Each  eigenvalue  has  finite  multiplicity  and  it  is  assumed  that  in 
the  enumeration  (5.4)  each  eigenvalue  is  repeated  according  to  its 
multiplicity.  There  exists  a corresponding  orthonormal  set  (4>.(x)} 
of  eigenfunctions  which  is  complete  in  L2  (G)  . Each  <}>.:  satisfies 


j e D(A(G) ) = L^(A,G)  and  A(G)<f)j  = Aj<J>j 
solutions  of  the  eigenvalue  problem 


<P 


Formally,  the  <£j  (x)  are 


^ ~ -f  + X 4>  = 0 for  x e G 

oX  j dX^ 

-1^-  = 0 for  x£  9G 


(5.6) 


Of  course,  if  3G  is  not  smooth  then  the  boundary  condition  is  the 
generalized  Neumann  condition  defined  in  section  4.  It  is  known 
that  the  first  eigenvalue  X0  = 0 is  simple  with  normalized 
eigenfunction 


const . 


(5.7) 


where  |g|  is  the  Lebesgue  measure  of  G. 


5.2  The  eigenfunction  expansion 

The  eigenfunctions  of  A may  be  constructed  by  separation  of 
variables.  From  a more  sophisticated  point  of  view,  A is  a sum  of 
tensor  products 


A = A(C)  € 1 + 1 8 A(R+) 


(5.b) 


where  R+  = (y:  y > 0}.  It  follows  that  the  eigenfunctions  of  A 

are  products  of  eigenfunctions  of  A(G)  and  A(R+) . The  spectral 
analysis  of  A(R+)  is  given  by  the  Fourier  cosine  transform  in 
L2(R+): 

' - rM 

cos  py  f (y)dy  (5.9) 


f (p) 

= L2(R+)-lim 

M-+CO 

f(y) 

= L2 (R+)-lim 

M-*00 

2l 

w 


1/2 

1/2 


cos  py  f(p)dp 


(5.10) 
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,f'W  ' '"MM 


(5.11) 


It  follows  that  a complete  normalized  family  of  generalized  eigen- 
functions for  A is  defined  by 


Wj(x,y,p)  = 


r 2^  1/2 


cos  py  <Mx),  p e R j = 0,1,2, ••• 

J (5.12) 


The  Plancherel  theory  for  A(G)  and  AC^),  quoted  above,  implies 
that 


f . (p)  = L2  (ELi.)  — lim 
J 


M-+°o 


r M 


0 J G 


Wj(x,y;p)  f (x,y)dxdy 


(5.13) 


exists  for  all  f G L2(S),  and  the  operator  , : L2(S)  -*■  L?  (R^.) 

defined  by  4>jf  = f^  has  range  <tjL2(S)  = L2(R+).  Moreover 

OO 

a f ii T2  ...  = y if.i.2  . 

L2(S)  l j l2(r+) 


(5.1/*) 


and 


j=0 


f (x,y)  = L2 (S) -lim  £ 


M,N-k» 


j=0 


M 


wj (x,y , p)  f j (p) dp 


(5.15) 


The  relations  (5.13)  and  (5.15)  are  frequently  written  in  the  more 
concise  symbolic  form 


£j(p)  ‘ 


wj  (x , y , p)  f (x ,y)dxdy 


(5.16) 


and 


f(x,y)  = l 


j=0  R+ 


Wj(x,y,p)  f j (p)dp 


(5.17) 


but  must  be  understood  in  the  sense  of  (5.13)  and  (5.15), 


Note  that,  formally,  f.(p)  is  just  the  L2 (S)  inner  product  of 
f(x,y)  and  the  eigenfunctiofl  (5.12).  For  a more  detailed  discussion 
of  this  expansion  see  [21]. 

The  generalized  eigenfunctions  (5.12)  are  locally  in  D(A)  and 
satisfy 
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AWj(*,*,p)  = -Aw_.  ( • , • ,p)  = (p2  + ^)Wj(*,*,p)  (5.18) 

This  fact  and  the  Plancherel  theory  imply  the  following  construc- 
tion of  the  spectral  family  of  A. 


5.3  Theorem 

If  {I1(X),X  > 0}  denotes  the  spectral  family  of  A = A(S)  then 
II(X)  has  the  eigenfunction  expansion 


Il<X)  f(x,y)  = I 


x.<x 

J- 


/x^XT 

J 


Wj(x,y,p)  f j (p) dp 


= l 


x.<x 

J- 


■ A-Aj 

2 

. 

0 

1/2 


cos  py  f j (p)dp 


(5.19) 


Vx) 


for  all  X > 0.  In  particular,  A is  an  absolutely  continuous  oper- 
ator whose  spectrum  is  the  interval  [X0  ,°°)  = [0,°°). 


Note  that  the  sum  in  (5.19)  is  actually  finite  by  (5.5). 
(5.19)  implies  that  the  eigenfunction  expansion  (5.17)  defines  a 
spectral  representation  for  A in  the  sense  of  the  following 
corollary . 


5.4  Corollary 

If  H'(X)  is  any  bounded  Lebesgue-measurable  function  of  X > 0 
then  for  all  f e L,,(S) 


N 


rM 


'i'(A)  f (x,y)  = L2 (S)-lim  £ 


3=0 


w (x,y,p)  'V  (p2  + X ) f (p)dp 

J J J 

(5.20) 


The  eigenfunction  expansion  (5.17)  defines  a decomposition  of 
the  Hilbert  space  L2(S).  To  describe  its  properties  let  f 6 L2 (S) 
and  define 


Pj  f(x,y) 


i.(x')  f(x',y)dx' 


<t>j(x) 


(5.21) 


f j (y)  4>j  (x) , j = 0,1 ,2, • • • 


where 


fj(y)  = | 4>j(x')  f(x',y)dx',  j = 0,1,2,' 


(5.22) 


The  orthonormality  of  {<}).}  in  L2(G)  implies  that  {P  : j •»  0,1 ,2,  • • •} 
defines  a complete  famil^  of  orthogonal  projections^in  LZ(S): 


p j = pj>  pjpk  - Vpk for  j’k  = °’1’2’ 


(5.23) 


l Pj  = 1 (3.24) 

j=0 

Moreover,  a simple  calculation  gives 

f(A)  P f ■ M(A)f  « w.  ( •,  •,p)'i'(p2  + X )f  . (p)dp  (5.25) 

J J R J J J 

for  j = 0,1,2, ••*.  In  particular, 


P.f(x,y)  = w (x,y,p)  f . (p) dp 


(5.26) 


An  equivalent  operator-theoretic  representation  is  P^  = 

JCj  = PjL2  (S)  = { f (x,y)  = fj(y)<t>j(x):  f G L2(R+)}  (5.27) 

then  (5.23)  - (5.25)  imply  the 


5.5  Corollary 

The  direct  sum  decomposition 


l2(s)  = l ex. 


(5.28) 


is  a reducing  decomposition  for  A. 

Note  that  each  3(.  is  isomorphic  to  L2(R  ) under  the  mapping 
f (x,y)  -*•  fj(y)  defined  by  (5.22). 


5.6  Solutions  in  L2(S)  of  the  propagation  prot lem 

Only  the  case  where  f 6 L2(S)  and  g € D(A-1^2)  will  be  dis- 
cussed. As  in  sections  3 and  4,  the  solution  in  L2(S)  of  the 
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propagation  problem  (4.1)  - (4.3)  has  the  form 
u(t,x,y)  = Re  {v(t,x,y)l 


(5.29) 


where 


v(t,*,“)  = exp  (-it  A1/2)h,  h = f + iA~l/2ge  L2(S)  (5.30) 

The  decomposition  (3.28)  implies  that 

OO 

v(t,x,y)  = £ Vj (t,x,y)  in  L2(S)  (5.31) 

j=0 

where 


Vj  (t ,x,y)  = PjV(t,x,y)  = v..  (t.y)^  (x)  6 L2(S)  (5.32) 

with 

(2)  1/2  f -itw.(p)  „ 

vj(t»y)  = hi  I cos  py  e J h^  (p)  dp  (5.33) 

R+ 

and 

u>.(p)  « (p2  + A.)  1/2  > X.l/2  > 0 (5.34) 

J J J 

In  the  theory  of  waveguides  (5.31)  is  called  a modal  decomposition 
and  the  partial  waves  v.=  (t,x,y)  are  called  waveguide  modes. 
Vj(t,x,y)  will  be  said  to  be  in  mode  j of  the  waveguide  S.  In 
particular,  mode  0 

v0(t,x,y)  = v0 (t,y)/ |G| 1/2  (5.35) 


will  be  called  the  fundamental  mode  of  S.  It  is  not  difficult  to 
show  that 


u0(t.y)  = Re  (v0(t,y)} 


= 7 {f0(y-t)  + *f0(y+t)>  + \ 


y+t 

g0 (y')d^ ' 

■ y-t 


(5.36) 


where  f0(-y)  = f0(y)  and  g0 (-y)  = g0 (y) . Note  that  the  modal 
waves  propagate  independently  in  the  sense  that,  different  modes 
are  orthogonal  in  L2(S)  for  all  t. 


The  spectral  representation  (5.31),  (5.32),  (5.33)  will  now 
be  used  to  study  the  asymptotic  behavior  for  t -*  00  of  solutions 
in  L2  (S ) . Because  of  the  independence  of  the  modes  it  will  be 
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enough  to  study  the  individual  modal  waves  (5.33).  The  substitu- 
tion 2 cos  py  = exp  (ipy)  + exp  (-ipy)  gives  the  decomposition 


Vj(t,y)  = v|(t,y)  + vT (t ,y) 


(5.37) 


where 


vl'(t.y)  = vT  (t  ,-y)  = 


(2tt)i/2 


i(yp-tu)  (p))  A 
e h^(p)dp 


(5.38) 


and  the  integral  converges  in  L2(R+)  (and  in  L2(R))  for  each 

h-  £ L2(R+).  The  special  case  of  the  fundamental  mode  is  discussed 

first. 


5.7  Asymptotic  wave  functions  for  the  fundamental  mode 

This  case  is  closely  related  to  that  of  section  3,  since 
w0(p)  = P f°r  all  P e Thus 


vo(t*y> 


ei(y-t)P  ho(p)dp  = Gfl(y-  t)  (5.39) 


where 


c„(y)  = 


(2tt)  1/2 


eiyP  hfl(p)dp  £ L2(R) 


(5. AO) 

Moreover,  it  is  easy  to  verify  by  direct  calculation  that  v0(t,y) 

(5.A1) 


= vj(t,-y)  = GQ(-y-t)  0 in  L2(R4.)  when  t ■*■  <*>.  Thus 


v0(t,y)  = G0 (y-t) 
is  an  asymptotic  wave  function  for  v0(t,y)  in  L2(R+): 

lim  lv  (t,*)  - v“(t,*)».  /„  = 0 (5.42) 

t-*» 

for  all  h„  £ L2(R+). 

For  the  higher  order  modes  j > 1 the  functions  ui.(p) 

* (p2  + A )V2  with  Aj  > 0.  For  these  cases  the  spectral  integrals 
(5.38)  all  have  the  same  form,  differing  only  in  the  value  of  Aj 
and  the  function  hj  £ L2(R+).  The  asymptotic  behavior  of  these 
integrals  may  be  determined  by  the  method  of  stationary  phase,  as 
follows . 
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5.8  Application  of  the  method  of  stationary  phase 
Consider  the  wave  function  defined  by 

v(t,y,A,h)  = (2tt)_  1/2  [ exp  {i(yp  - tu)(p,  A))  }h(p)dp  (5.43) 

^ R+ 
where 

w(p, A)  = (p2+  A)  1/2  > A 1/2  > 0 

and 

h G L2(R+) 

The  phase  function 

©(p.A.y.t)  = yp  - tw(p, A) 
is  stationary  with  respect  to  p if  and  only  if 
30(p,A,y,t)/3p  = y - t3w(p, A)/3p  = 0 
or 

t 3p  “ U(p,A;  (d2+A) 'ft 


(5.44) 


(5.45) 


(5.46) 


(5.47) 


(5.48) 


The  function  U(p,A)  defined  by  (5.48)  is  the  group  velocity  [5] 
for  the  wave  function  (5.43).  Note  that 


3U (p, A)  _ 32lo(p,A) 


> 0 


(5.49) 


3p  3p2  (p?+A)3/2 

and  hence  U(p,A)  is  a monotone  increasing  function  of  p.  Moreover, 
0 < U(p,A)  < 1 for  all  p > 0 and  A > 0 (5.50) 

Hence  for  t > 0 equation  (5.48)  has  the  unique  solution 
iR/jLl2 


P = 


[l-  (y/t); 


U/2 

y2*.. 

[t2-  y2 

1/2 


> o 


if 


0 < y/t  < 1 


(5.51) 


(5.52) 


and  has  no  solution  for  other  positive  values  of  t.  The  principle 
of  stationary  phase  asserts  that  for  large  values  of  y‘ + t2  the 
stationary  point  (5.51)  will  make  a contribution 


w 
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v°°(t,y,  A,h) 


0i(yp-tw(p,A)-ir/4) 

-e- h(p),  (5.53) 

(t3U(p,A)/3p) 1/2 


to  the  integral  (5.43),  where  is  the  characteristic  function 
of  the  set  (5.52).  More  precisely,  if  h € P(R^)  then  the  following 
error  estimate  is  known  [2,23]. 


5.9  Theorem 

Let  h £ V(R+)  and  define  the  remainder  q(t,y,A,h)  by 
v(t,y,A,h)  = v (t,y,A,h)  + q(t,y,A,h)  (5.54) 

Then  there  exists  a constant  C = C(A,h)  such  that 

|q(t,y,A,h) | < C/ (y2  + t2)  3^*  for  all  y £ R and  t > 0 (5.55) 

It  follows  from  (5.54)  and  (5.55),  by  direct  integration,  that 
for  all  h £ P(R.) 

T 

lim  llv(t,*,A,h)  - vC°(t,  • , A ,h)ll  , . = 0 (5.56) 

t-K» 

The  stationary  phase  method  is  not  applicable  to  (5.43)  when  A = 0. 
However,  the  results  for  this  case  are  described  by  the  same 
equations  if 


v°°(t,y,0,h)  = ( 2tt  ) — 1/2 


exp  (iyp)  h(p) dp 


(5.57) 


With  this  notation,  (5.56)  with  A = 0 is  equivalent  to  (5.42). 


The  estimate  (5.55)  implies  (5.56)  for  all  h £ P(R , ) . For 
more  general  h £ L2(R+)  the  estimate  (5.55)  may  not  hole.  Never- 
theless, the  following  results  hold. 


5.10  Theorem 

For  all  A > 0 and  all  h £ L2 (R+) 
v°°(t,*,A,h)  £ L2  (R+)  for  all  t ^ 0 


(5.58) 


I 

I 


t -*■  v (t,*,A,h)  6 L2  (R+)  is  continuous  for  all  t ^ 0 (5.59) 


H v (t,*,X,h)llL  ^ < II h II L ^ for  all  t ^ 0 


(5.60) 


Moreover,  the  relation  (5.56)  holds  for  all  h E L2(R+). 

Properties  (5.58)  - (5.60)  follow  from  the  definitions  (5.53), 
(5.57)  by  direct  integration.  Moreover,  the  validity  of  (5.56)  for 
all  h € L2 (R . ) follows  from  the  special  case  h 6 ^(R^) , the  density 
of  P(R_|_)  in  L2 (R+)  and  the  uniform  boundedness  in  t of 
Kv(t,*,A,h)ll  . . and  II v (t,*,A,h)ll1  , ..  More  detailed  proofs 

may  be  found  in  [22,40]. 


5.11  Asymptotic  wave  functions  for  the  higher  order  modes 
Define  the  modal  asymptotic  wave  functions  by 
v“(t,y)  = v°°(t,y,Aj,hj) , j = 0,1,2,*** 

Then  (5.38),  (5.43)  and  (5.56)  imply 


lim  II  v (t,  •)  - v (t,  *)H 


r ' ' l2(r+) 


= 0,  j = 0,1,2, 


(5.61) 


(5.62) 


Moreover,  (5.38)  for  v^  and  (5.43)  imply 
(t,y)  = v(t,-y,A^,hj) 


= (2tt)  — 1/2  exp  {~i(yp + tw(p,A  )}  h (p)dp 

J T)  J J 


(5.63) 


The  stationary  phase  method,  applied  to  (5.63),  implies  that 


lim  Bv  (t,-)||  , . = 0 

t-+oo  J l2(r+) 


(5.64) 


because  the  phase  yp  + ta)(p,A.)  in  (5.63)  has  no  stationary  points 
when  y > 0 and  t > 0.  Combining  (5.37),  (5.62)  and  (o.64)  gives 


lim  «v  (t,‘)  - v"(t,*)ll,  \ = 0 ( 

t-MX>  J J 

for  all  h.  € L2(R  ) and  j = 0,1,2, ••*.  The  results  and  the 
decomposition  (5.31),  (5.32)  imply  the 


(5.65) 


J 


Then 


and 
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5.12  Asymptotic  convergence  theorem 
For  all  h e L2(S)  define 


v°°(t,x,y)  = l v“(t,y)  4)^  (x) , (x,y)  e S (5.66) 

j=0 

v(t,*,*)  e L2(S)  for  all  t * 0 (5.67) 

t -*■  v(t,  *,  •)  £ L2(S)  is  continuous  for  all  t ^ 0 (5.68) 

0 v°°(t,  • , •)!  L ^ 5 U hll  L ^ for  all  t ^ 0 (5.69) 


lim  llv(t,*,*)  - v (t,«,-)ll  , . = 0 (5.70) 

t-«o  L2'*) 

The  proof  of  this  result  will  be  outlined.  First,  note  that 
the  convergence  in  L2(S)  of  the  series  in  (5.66)  follows  from  the 
orthogonality  of  its  terms  in  L?(S),  (5.60)  which  implies 

,Vj(t,,)VL2(S)  = * vj  ( t > * > !l2  (R_^)  - *Vl2(R+)  (5,71) 

for  all  t ^ 0 and  (see  (5.14)) 


,hlMs>  ' I "Vmv 

j-o 


< 00 


(5.72) 


Properties  (5.68)  and  (5.69)  follow  from  (5.59)  and  (5.60),  applied 
to  Vj(t,y).  Finally,  to  verify  (5.70)  note  that  for  j = 0,1,2, 

,Vj(t,  *>  - Vj(t,  *)I*L2 (R+)  - ‘)!IL2(R+)  + l,Vj(t,‘)llL2(R+) 


(5.73) 


< 20  h 


J MR+) 


for  all  t 0.  It  follows  that 


lv(t,*,*)  - v (t, •,•)!!* 


= l »Vj(t,*)  - vj(t.*)|,L2(R  ) 

j=0 


(5.74) 


44 


N *>  (5.74  Cont) 

< l lVj(t,-)-v”(t,-)l^(R  } + 4 l ) 

j=0  j =N+1 

for  N = 0,1,2,  •••.  Fixing  N and  making  t -*•  00  gives,  by  (5.65) 

00 

lim  lv(t,-,0  - v“(t,*,-)l*a(s)  £ 4 l ) (5.75) 

t’>CO  ~=N+1  + 


for  N = 0,1,2, ••*.  Thus  (5.70)  follows  from  (5.72)  and  (5.75). 

If  f £ L2(S)  = D(A!/2)  and  g £ L2(S)  then  the  same  method 
can  be  used  to  show  convergence  in  energy: 


w 

lim  E(u  - u ,S, t)  = 0 

£-X» 


(5.76) 


where  u (t,x,y)  = Re  {v  (t,x,y)}  but  the  details  will  not  be 
recorded  here. 


6.  SCATTERING  BY  OBSTACLES  AND  JUNCTIONS  IN  TUBULAR  WAVEGUIDES 

The  analysis  of  section  5 is  extended  to  compound  tubular 
waveguides  in  this  section.  The  mathematical  problem  is  the 
initial-boundary  value  problem  (4.1)  - (4.3)  for  an  unbounded 
domain  !1  C R3  of  the  form 


ft  = ft.  U S U ...  U s 

0 l m 


(6.1) 


where  ft0  is  a bounded  domain  and  S1,,**,S  are  disjoint  uniform 
semi- inf inite  cylinders.  Examples  include  waveguides  of  the  types 
described  at  the  beginning  of  section  5 and  many  others.  It  will 
be  assumed  that  ft  £ LC. 


6.1  Notation 


It  will  be  convenient  to  think  of  R3  as  a 3-dimens.ional 
differentiable  manifold.  The  generic  point  of  R3  will  be  denoted 
by  q.  A special  Cartesian  coordinate  system 


/CL  CL  CL v - / u vav  D 

(Xj ,x2,y  ) - (x  ,y  ) £ R 


a ax 


(6.2) 


may  be  associated  with  each  semi-infinite  cylinder  (a=l,**",m) 
in  such  a way  that 

Sa  = {q  £ R3:  xa(q)  £ G^  and  ya(q)  > 0}  (6.3) 


(6.3) 
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where  G C R is  a bounded  domain.  The  assumption  that  Q € LC 
implies  that  Ga  6 LC  for  cx  = l,***fm  and  hence  that  each  A(Ga)  has 
a discrete  spectrum  with  eigenvalues 


A „ < A , 
ao  - ai 


< A 


OU 


such  that 


(6.4) 


lim  A „ = 00  (6.5) 

£-*»  a* 

and  corresponding  eigenfunctions 

*a0(x“)  = 1/IgJ1'2,  4>ai(xa),  <t>a2(xX)  (6.6) 

which  form  a complete  orthonormal  sequence  in  L2 (G^) . 


6.2  Solutions  of  Aw  = Aw  in  Sa 

Suppose  that  w is  locally  in  D(A);  i.e.,  4*w  £ D(A)  for  every 
4>  e D(R2).  Then  the  completeness  of  the  eigenfunctions  (6.6) 
implies  that 

00 

w(q)  = I wa£(y0l)  *as,C*a>  for  311  Q e sa  <6-7> 

£=0 

where  xa  = xQ(q),  ya  = ya(q) • Moreover,  if 

Aw  = Aw  in  S^  (6.8) 

a 

then  the  coefficients  wa£(ya)  will  satisfy 

wa<t(ytX)  + (X  " wct£(y01)  = 0 for  a11  yCl  > ° (6’9) 

In  particular,  if  it  is  assumed  that 

A ^ A a = 1,  • • • ,m;  i = 0,1,2, •••  (6.10) 

then 


V(y0l)  = Ca £ exp  {i/XrAaJl  y0t}  + CaA.  exp  {_i/X~Xa!l  yCt}(6*11> 

where,  for  definiteness,  p1^2  > 0 for  p > 0 and 


►OTTA 


a£ 


f ‘x  - V £"  X > V 

for  A < 


(6.12) 
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6.3  Eigenfunctions  of  A and  non-propagating  modes 


It  was  discovered  by  F.  Rellich  [27]  that  the  operators  A for 
waveguide  regions  of  the  form  (6.1)  may  have  a point  spectrum.  A 
point  X G R is  in  the  point  spectrum  of  A if  and  only  if  there  is 
a non-zero  function  w G D(A)  such  that  Aw  = Xw.  In  particular, 
the  requirement  that  w e L2(ft)  implies  that  in  the  expansions 
(6.7),  (6.11)  the  coefficients  = C~^  = 0 for  X > Xa£  and 
= 0 for  X < X^£.  Thus  any  eigenfunction  of  A must  have  the 

form 


w(q)  = l c+£  exp  {-(Xai- X)l/2  y°}  ^(x3)  (6.13) 


ai 


aH 


U:X<Aae} 


for  all  q G S . In  particular,  the  eigenfunctions  are  exponentially 
damped  in  each  cylinder  S . 


D.  S.  Jones  [17]  has  shown  that  the  point  spectrum  of  A is  a 
discrete  subset  of  ( 0 , °°) ; i.e.,  each  eigenvalue  has  finite  multi- 
plicity and  each  finite  subinterval  of  (0,°°)  contains  at  most  a 
finite  number  of  eigenvalues.  Thus  if  the  point  spectrum  of  A is 
not  empty  then  there  exists  an  M such  that  1 < M < 00  and  ^(n) > 

1 < n < M,  is  an  enumeration  of  the  eigenvalues  of  A,  each  repeated 
according  to  its  multiplicity.  It  may  be  assumed  that 


0 < \n)  - X(n+1)  for  1 - n < n+l  * M (6.14) 

The  corresponding  eigenfunctions  will  be  denoted  by  w,  . . The 
subspace  spanned  by  {w^  : 1 < n < M}  will  be  denoted  by  ^(A) 
and  called  the  subspace  of  discontinuity  of  A [18].  Thus 


i 


3^(A)  = {w  = l en  w (n):  l | cn  | 2 < °°}  (6.15) 

l<n<M  l<n<M 

It  is  known  that 

L2(ft)  = 'K^(A)  © K0 (A)  (6.16) 

where  ^(A),  the  orthogonal  complement  of  7^(A)  in  L2(ft),  is  that 
largest  subspace  of  L2(ft)  on  which  the  spectral  measure  of  A is 
continuous.  ^(A)  is  called  the  subspace  of  continuity  of  A [18]  . 
Moreover,  (6.16)  is  a reducing  decomposition  for  A [18]. 

If  the  initial  state  of  an  acoustic  field  in  ft  satisfies 
u(0,‘)  = f G 3(P(A)  and  9u(0,*)/9t  = g G 3(P(A)  then  h = f + i.A~1^2g 
G J{P(A)  and  hence 
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v(t,q)  = exp  (-itA1'2)  h(q) 


l c(n)  exp  (-itX&)  w(n)(q) 
l<n<M 


It  follows  that  the  energy  of  the  acoustic  field  u( t , q)  = Re{v(t , q) } 
in  any  bounded  portion  of  0.  is  an  oscillatory  function  of  t.  In 
particular,  there  is  no  propagation  of  energy  in  the  cylinders  S^. 
For  this  reason  the  eigenfunctions  W(n) (q)  are  called  non- 
propagating modes  of  the  waveguide.  By  contrast,  it  is  shown 
below  that  for  fields  with  initial  state  in  3£"(A)  the  energy  in 
every  bounded  portion  of  f2  tends  to  zero  when  t ■*  «>  and  hence  all 
the  energy  propagates  outward  in  the  cylinders  S^. 


6.4  Generalized  eigenfunctions  of  A 

The  operator  A has  two  families  of  generalized  eigenfunctions, 
analogous  to  the  functions  w+(x,p)  and  w_(x,p)  of  section  4,  each 
of  which  spans  the  subspace  ^(A) . The  structure  and  properties 
of  these  functions  are  described  next. 

Consider  a single  term  in  the  expansion  (6.7)  for  the  cylinder 
Sa.  It  has  the  form  (cf.  (6.11)) 

wa£(q)  = Ki exp  yCX} 

(6.18) 

+ c;£  exp  ya})  *ae(xa) 

Cl  QJ 

where  q **  (x  ,y  ).  Assume  that  X > X £,  so  that  (6.18)  represents 
a propagating  mode  in  S^,  and  write 

P = (X  - Xa^2  > 0 (6.19) 

and 

A'«  = V(p)  - (p2  + A^)1'2  > A •«  (6.20) 

If  one  associates  a time-dependence  exp  { — iX 1 2t}  = exp  {-iwai,  (p)  t} 
with  (6.18),  as  in  the  spectral  representation  of  v(t,*) 

= exp  (~itA^2)h,  then 

wQjl(q)  exp  *_lu,a£(p)t*  = Ca£exp  ^i(pya~  wa£(p)t)^af  (xa) 

(6.21) 

+ Ca£  exp  ^-i(I,ya+u)a£(p)t)^a.V.(xa) 
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is  the  sum  of  an  outgoing  wave  in  Sa,  with  coefficient  CQ£,  and  an 
incoming  wave  with  coefficient  C^£.  For  this  reason,  a solution 
of  (6.8)  of  the  form 

CaJl  exp  (ipy0l)  (6.22) 

will  be  called  an  "outgoing"  wave  in  in  mode  Z,  while  a solution 
of  the  form 

CaJl  GXp  (_iPya)  ^(x*)  (6.23) 

will  be  called  an  "incoming"  wave  in  Sa  in  mode  Z.  Note  that  this 
terminology  is  based  on  the  convention  that  the  time-dependence  is 
exp  (-iwa^(p)  t) , as  in  (6.21).  If  a time-dependence  exp  (icoap(p)t) 
were  used  it  would  be  necessary  to  interchange  the  terms  "outgoing" 
and  "incoming." 

In  the  case  of  the  uniform  semi-infinite  cylinder  of  section 
5,  m is  equal  to  1 and  the  generalized  eigenfunctions  have  the  form 


w^U.y.p) 


1 

(2 it)  »/* 


exp  (ipy)  <P^(x) 


+ (2tt)V2 


exp  (-ipy)  4>^(x) 


(6.24) 


Thus  they  are  the  sum  of  an  incoming  and  an  outgoing  wave  in  mode 
Z,  with  equal  amplitudes  and  phases.  This  symmetry  is  due  to  the 
symmetry  of  the  waveguide.  In  the  general  case  of  a compound 
waveguide  (6.1)  it  is  possible  to  prescribe  the  amplitudes  and 
phases  of  the  incoming  (resp.,  outgoing)  waves  in  each  cylinder  Sp 
and  mode  Z.  The  amplitudes  and  phases  of  the  outgoing  (resp., 
incoming)  waves  in  each  cylinder  Sg  and  mode  m are  thereby  deter- 
mined. The  most  useful  generalized  eigenfunctions  are  those  that 
have  an  incoming  (resp.,  outgoing)  wave  of  prescribed  amplitude 
and  phase  in  a single  prescribed  cylinder  S and  mode  l.  They  may 
be  described  as  follows. 


6.5  Definition 

The  mode  (a, £)-outgoing  eigenfunction  for  ft  is  the  function 
w^(q,p)  defined  by  the  properties 

w^(-,p)  is  locally  in  D(A)  (6.25) 

(A  - w^(p))  w^(q.p)  - ”(A  + wa£(p))wajl(q,p)  = 0 (6.26) 


for  all  q € (1  and 
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w«8.(q’p)  = '(2~Tffi/2  exP  (-ipy°)  4>„j(xa) 


(6.27) 


+ l Ca*,Bm(p)  6XP  {i/p2  + Xa£-A3m 


m=0 


for  all  q 6 Sg  (B  = l,2,*‘,,m).  Similarly,  the  mode  (a, £)-incoming 
eigenfunction  for  is  the  function  wZo(9>p)  defined  by  the  prop- 
erties that  w~£(*,p)  is  locally  in  D(A),  (A+ u)2£(p)  )w~  (q,p)  = 0 
for  all  q£  !1  and,  for  q e (B  = 1 , 2,  • • • ,m) 


wa*(q’p)  =T2^-exp  (ipya)  4>aJl(xa) 


(6.28) 


+ 2 C^,Bm(p)  6XP  {-i,/p'  + Xa£-  XB«n 


m=0 


where  _i/p2 + Xa)T  ABm  < 0 f0r  XBm  > Xa*  + p* ' 

The  eigenfunction  wj^(q,p)  may  be  interpreted  physically  as 
the  steady-state  acoustic  field  in  the  waveguide  (2  due  to  a single 
incoming  wave  (6.23)  in  cylinder  S and  mode  £,  with  amplitude  and 
phase  defined  by  C^(p)  = l/(2ir)1'  , and  no  incoming  waves,  in  the 
other  cylinders  or  in  the  other  modes  of  cylinder  S^.  The  ampli- 
tudes and  phases  of  the  corresponding  outgoing  waves  are  defined  by 
the  coefficients  gm(p)  which  are  determined  by  the  incident  wave 
and  the  geometry  of  Note  that,  in  general,  an  incoming  wave  in 
mode  (a,£)  will  produce  outgoing  waves  in  all  the  cylinders  and 
modes;  i.e.,  scattering  produces  coupling  among  the  cylinders  and 
modes . 

g 

The  form  of  the  exponential  which  multiplies  4>gm(xP)  in  (6.27) 
is  determined  by  the  requirement  (6.26).  Note  that  the  sum  in 
(6.27)  includes  propagating  modes  with  Agm  < Xa£  + p‘  and  modes 
"beyond  cutoff"  with  Agm  > Xa£  + p2.  The  latter  decrease  exponen- 
tially when  y®  -*■  00 . 


The  eigenfunctions  w^(q,p)  have  an  interpretation  analogous 
to  that  of  wjj^(q,p),  but  with  "outgoing"  and  "incoming"  inter- 
changed. It  is  easy  to  verify  from  the  defining  conditions  that 
the  two  families  satisfy  the  relation 


wa£(q,p)  = wa£(q,p) 


(6.29) 


M 


t 


t 


3 
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The  case  of  the  uniform  semi-infinite  cylinder  is  a very  spe- 
cial case  in  which  m = 1 (so  that  no  index  a is  needed)  and 

w^(q,p)  = w“(q,p)  = w£(x,y,p),  q ~ (x,y)  (6.30) 

(see  (6.24)).  Moreover,  in  this  case  the  symmetry  implies  that 
there  is  no  coupling  between  different  modes: 


Existence  and  uniqueness  theorems  for  the  eigenfunctions 
w^£(q,p)  were  proved  in  [21].  The  following  notation  will  be  used 
to  formulate  them. 


+ Ve  VA)> 

where  Op(A)  = {A  = A^^:  1 < n < M} . Similarly, 

W!'£V  P*  + °<A<G(j»> 

where  o(A(G^))  = {X  = X^:  £ = Finally 

ID 

Zai  = zal^  U U zaf,(Gg) 

0=1 

and 

m 00 

z - 1)  IK*. 

a=l  2=1 


(6.32) 


(6.33) 


(6.34) 


(6.35) 


Note  that  the  information  on  the  spectra  of  A(Ga)  and  A given 
above  implies  that  each  of  these  sets  is  a denumerable  subset  of 
R+.  The  results  of  [21]  imply  the  following  theorem. 


6.6  Theorem 


Let  0 € LC  be  a waveguide  domain  of  the  form  (6.1). 
each  p e R+  - Z,  each  a = l,***,m  and  each  2 = 1,2, •••  the 
functions  W^£(*«P)  and  exist  and  are  unique. 


Then  for 
eigen- 


6.7  The  eigenfunction  expansion  theorem 

The  families  KjK’.p):  p C-  R+-Z;  a=l,**,,m;  2=0,1, 2, •••} 
and  {w~£(*,p):  Pe  R+-Z;  a = !,••*, m;  2 = 0,1 ,2,  • • •}  define  two 
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complete  sets  of  generalized  eigenfunctions  for  the  part  of  A in 
the  subspace  of  continuity  3^ (A) . The  eigenfunction  expansions, 
which  are  of  the  Plancherel  type  described  in  the  preceding  sec- 
tions, may  be  formulated  as  follows. 

6.8  Theorem 
Define 

Sa,M  = {q  6 r3;  *a(q>  e Ga  and  0 < yCX(q)  < M}  (6.36) 


and 


ft.  = (2.  us,  M u •••  u s M 

M 0 l ,M  m,M 

Then  for  all  f 6 L2(f2)  the  limits 


(6.37) 


i 


faJl(p)  = L2  (R+) -lim 


M-+00 


wajl(q,p)  f (q)  dV 


(6.38) 


(2. 


M 


exist,  where  dV  is  the  element  of  Lebesgue  measure  in  R3 . ^More- 
over, the  operators  L2(fi)  -*■  l^CR^.)  defined  by 

have  range  L2 (R+)  and,  if  Pc  denotes  the  orthogonal  projection  of 
L2  (f2)  onto  -K0 (A)  then 


m 


I Pcf B 2 


L2(S)  lfa£"L2(i4) 

a=l  £.=0 


(6.39) 


for  all  f e L,(£2),  and 


m N 


PCf(q)  = L2(f2)-lim  l l 
M’N~  a=l  £=0 


M 


wa^(q»p)  frtJ(p)dp  (6.40) 


a£ 


The  relations  (6.38)  and  (6.40)  will  be  written  in  the 
symbolic  form 


K £(p) 


aSL 

(2 

m 00 

Pcf(q)  = I l 
a=l  £=0 


v’  n(q,p)  f(q)  dv 


+ ± 
wa£(q,p)  f«0(p)dp 


a£ 


(6.41) 


(6.42) 


■4 


but  they  must  be  understood  in  the  sense  of  (6.38),  (6.40). 
following  corollaries  are  almost  immediate;  see  [13,21]. 


6.9  Corollary 


For  each  f £ L2(f2)  the  limits 


f*£(q)  = Mn)-lijn  wa£(q,p)  ^a£(p)dp 


exist  and 


(fa£,f8m)L2(S2)  = ° whenever  <a>£>  * 


Moreover, 


m 00 


pCf  = l l 4 


a=l  £=1 


6.10  Corollary 
Define 

3T£  = {f^  e L2(fi);  f€L?(Q)}  (6.4 

Then  each  is  a closed  subspace  of  K0  (A) , and  3f^m  are 

orthogonal  whenever  (a,£)  / (8>n>)  and 

m 00  m 00 

I'm  - l !®<r  ! I <6-' 

a=l  £=1  a=l  £=1 

The  eigenfunction  expansions  (6.40)  provide  the  following 
construction  of  the  spectral  family  of  A in  3Cc(A). 


6.11  Theorem 

If  {11(A):  A > 0}  denotes  the  spectral  family  of  A then 


II(A)Pcf(q)  = i l 


wa£(q’p)  ?a£(p)dp 


a=l  Xa£<X  1 o 


53 


for  all  A 2 0.  In  particular,  APC  is  an  absolutely  continuous 
operator  whose  spectrum  is  [O,00). 

Note  that  the  sums  in  (6.48)  are  actually  finite  because 
Xqji  -*■  00  when  Z -*■  °°.  (6.48)  implies  that  the  eigenfunction  expan- 

sions (6.40)  define  spectral  representations  for  A in  the  sense 
of  the  following  corollary. 


6.12  Corollary 

If  ’i'(A)  is  any  bounded  Lebesgue-measurable  function  of  X > 0 
then  for  all  f e ?(^(A)  = PcL2(ft) 

m 00 

4'(A)f(q)  = l l wa£(q,p)  T(pZ  + Xa£^a£(p)dp  (6,49) 

a=l  £=1  R+ 

It  follows  from  (6.47)  and  (6.49)  that  the  eigenfunction  ex- 
pansions (6.40)  define  reducing  decompositions  of  ^(A).  More 
precisely,  the  following  generalization  of  the  results  of  section 
5 is  valid  [13,21]  . 


6.13  Corollary 

+ + + 

The  operator  Pa£  defined  by  P~^f  = f“^  is  an  orthogonal  pro- 
jections of  L2(f0  onto  ^ and 

m °° 

pC  = l l pa£  (6'50) 

oi=l  1=1 


Moreover, 

P“^II(A)  = IT(A)P~^  for  all  A > 0 (6.51) 

and  hence  (6.47)  defines  reducing  decompositions  for  APC. 

+ The  surjectivity  of  4>“^:  L2(fi)  -+•  L2  (R  ) , the  definition  of 

P~^  and  (6.43)  imply  that  for  all  a = l,***,m  and  Z = 1,2, ••• 


(6.52) 


In  particular  the  eigenfunction  mappings  <t>a£  are  partial  isometries 
[18]  with  initial  sets  and  final  sets  L2 (R+) . 
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1 


6.14 


Solution  in  ^(A) 


of  the  propagation  problem 


Only  the  case  where  f 6 L2(ft)  and  g £ D(A-1^2)  will  be  dis- 
cussed here.  For  such  initial  states  it  follows,  just  as  in 
sections  3,  4 and  5,  that  the  solution  in  L2(ft)  of  (4.1)  - (4.3) 
has  the  form 


u(t,q)  = Re  {v( t , q) } (6.53) 

where 

v(t,*)  = exp  (-itA(ft)  ^2)h,  h = f + iA(ft)_1/2g  e L2(ft)  (6.54) 

Moreover,  the  case  where  h £ 3(^(A)  was  discussed  above.  Hence, 
only  the  case  where  h £ 3(C(A)  remains  to  be  analyzed.  In  this 
case  v(t,‘)  £ 3(C(A)  for  all  t £ R and  (6.47)  implies  that  v(t,q) 
has  the  decompositions 

m 00 

v(t,q)  = l l in  ^(A)  (6.55) 

a=l  £=1 


where 


va£(t’q) 


exp 


(-itA(fi)1/2)h^ 

+ 

wctf/q,p)  exp  (-ita)a£(p))  ha£(p)dp 


(6.56) 


and  wa£,(p)  is  given  by  (6.20).  The  two  decompositions  defined  by 
(6.55),  (6.56)  will  be  called  modal  decompositions,  in  analogy 
with  the  simple  case  of  section  5,  and  the  partial  wave  v^^(t,q) 
will  be  said  to  be  in  mode  (±,a,£)  of  the  compound  waveguide  ft . 
Note  that  for  the  uniform  semi- inf inite  cylinder  of  section  5 the 
(+,£.)  and  (-,£)  modes  coincide  (see  (6.30)). 


6.15  Transiency  of  waves  in  7^  (A) 

The  absolute  continuity  of  the  operator  A in  the  subspace 
^(A)  implies  that  all  waves  in  .^(A)  are  transient  in  the  sense 
of  the  following  theorem  [38,42]. 


6.16  Theorem 

If  ft  £ LC  is  a waveguide  domain  (6.1)  then  for  every  h £ (A) 

and  every  compact  set  K C R3 
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lim  0 exp 

£-*0O 


<-ltA'«)hls{M)  - 0 


(6.57) 


Thus  the  decomposition  L?(ft)  = 3{^(A)  ©J^A)  splits  every 
h G L2 (0)  into  a sum  of  a non-propagating  and  a propagating  state. 
In  particular  the  partial  waves 


Va£^t’’^  = exp  (-itAl/2)  p^h  6 3^ (A)  (6.58) 

and  hence  (6.57)  with  K = implies 

lim  »v*A(tf-)lL2(0#)  = 0 (6.59) 

for  a = 1,2, and  9,  = 0,1,2, •••.  Thus  waves  in  3^ (A)  ulti- 
mately propagate  into  the  cylinders  Sa.  The  eigenfunction  expan- 
sion for  A will  now  be  used  to  calculate  the  asymptotic  form  of 
these  waves. 


6.17  Asymptotic  wave  functions 

Let  h € J^CA)  and  consider  the  representation 


m 


v( t , * ) — exp  (-itA1/2)h  = £ £ v'^t,*) 

a=l  9,-0 


(6.60) 


defined  by  the  incoming  eigenfunctions  W^£(9»P)«  Substituting  the 
development  (6.28)  for  w~£(q,p)  in  Sg  into  the  integral  (6.56)  for 
v~^(t,q)  gives  the  representation 


v'c(t,q)  = 6 ov(t,ya,A  „ ,h  .)<J>  „(xa)  + v * (t.q). 


(6.61) 


q G Sf 


where  v(t,y,A,h)  is  defined  by  (5.43)  and 


v^U.q)  = l V^gm(t,y  ) ^6ra(x  ),  q e S 


a 9,,  3m 


3m 


(6.62) 


m=0 


with 
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Va*,em(t’y) 


This  equation  defines  for  the  case  where  Xgm  > Xa£.  In  the 

case  where  Xgm  < Xa£  the  first  integral  is  absent  and  the  second 
has  the  lower  limit  zero.  The  method  of  stationary  phase  may  be 
applied  to  show  that  v'p  p (t,y)  satisfies  an  estimate  of  the  form 
(5.55)  for  y > 0,  t > 0 because  the  integrals  in  (6.63)  have  no 
points  of  stationary  phase  in  this  region.  It  follows  that 

lim  B v^£(t, *)IL  ^ ^ = 0,  $ = 1, 2, • • • ,m  (6.64) 

The  proof  may  be  based  on  a convergence  lemma  like  that  of  section 
3.  The  details  will  not  be  given  here.  See  [22]  for  a more 
complete  discussion. 

In  (6.61)  the  term  v(t ,ya, X^^jh^g,)  has  the  form  (5.43)  studied 
in  section  5.  Hence,  if  v°°(t,y,X,h)  is  defined  as  iri  that  section 
and 


^Bm  }'al 


(-A*™  - Ko  -V 


exp 


* exP  (~itU)a£(P))  Ca£,3m(P)ha£(P)dp 

„ (6-63) 

e*P  + *6m  y+tajai(p))} 

^Bm  ^ exit 

* Ca£,Bm(p)  ^Jl(P)dP 


va*(t’q)  = v“(t’ya’\*r*W  q G Sa  (6-65) 

then  (6.64)  and  the  results  of  section  5 imply 


11« 


* ^afiva2(t ’ L2 (Sg) 


= 0, 

e - i,2, 


,m 


(6.66) 


In  particular,  + 0 in  L2(Sg)  for  t °°  and  all  6 t a; 

i.e.,  ^£(1,*)  is  asymptotically  concentrated  entirely  in  S^. 

The  asymptotic  wave  functions  for  v(t,q)  will  be  defined 


by 
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m 00 

v°°(t,q)  = I l X^q)  v^U.q).  qen  (6.67) 

a=l  A=0 

where  Xa^q)  denotes  the  characteristic  function  of  Sa.  Note  that 
the  terras  in  this  sura  are  orthogonal  in  L2(D)  by  (6.65).  The  de- 
compositions (6.60),  (6.67)  and  the  convergence  results  (6.59) 
and  (6.66)  imply  the 


6.18  Theorem 

If  SI  6 LC  is  a waveguide  domain  (6.1)  then  for  every  h e jfC(A) 
the  wave  function  v(t,»)  = exp  (~itA1'2)h  satisfies 


lim  I v( t , 

fc-xx> 


) - v (t,-)ll 


L,(fl) 


= 0 


(6.68) 


The  proof  is  essentially  the  same  as  for  the  special  case 
described  in  section  5.  The  convergence  in  energy,  when  h has 
finite  energy,  can  be  proved  by  the  same  methods. 


7.  PROPAGATION  IN  PLANE  STRATIFIED  FLUIDS 

The  propagation  of  localized  acoustic  waves  in  a plane  strati- 
fied fluid  which  fills  a half-space  is  analyzed  in  this  section. 

The  asymptotic  wave  functions  for  such  media  are  shown  to  be  the 
sura  of  an  asymptotic  free  (hemispherical)  wave  and  an  asymptotic 
guided  wave  which  propagates  parallel  to  the  boundary.  This 
structure,  which  is  intermediate  between  that  of  a homogeneous 
fluid  and  that  of  a tubular  waveguide,  is  called  an  open  waveguide 
in  the  physical  literature. 


7.1  Plane  stratified  fluids 

An  inhomogeneous  fluid  will  be  said  to  be  plane  stratified  if 
the  local  sound  speed  c(x)  and  density  p(x)  are  functions  of  a 
single  Cartesian  coordinate.  This  condition  can  be  written 

c(xj,x2,x3)  = c (x  3 ) | 

> (7.1) 

P(X1,X2,X3)  = p ( X 3 ) J 

with  a suitable  numbering  of  the  coordinates.  It  will  be  convenient 
to  denote  x3  by  a single  letter  and  write 
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x *=  (Xj,x2)  G R2,  y = x3  e R,  (x,y)  G R3 


(7.2) 


This  notation  is  used  in  the  remainder  of  this  section. 


7.2  Propagation  in  a stratified  fluid  with  a free  surface 
A stratified  fluid  filling  a half-space 


R+  = {(x,y):  x G R2  and  y > 0} 


(7.3) 


is  often  used  as  a model  in  the  study  of  acoustic  wave  propagation 
in  oceans  and  deep  lakes.  If  the  surface  {(x,0):  x € R2}  '3  free 
the  corresponding  initial-boundary  value  problem  is  (see  section  2) 


92u  2/  \ f <,  ^ 1 3u 

jp-  - c (y)p(y)  ST  ^ 

J ( J 


u(t,x,0)  = 0 for  t > 0,  x G R2 


= 0 for  t > 0, 


(x,y)  G R? 


(7.4) 


(7.5) 


u(0,x,y)  = f (x,y)  and  9u(0,x,y)/3t  = g(x,y) 


(7.6) 


for  (x,y)  G R| 

where  in  (7.4)  j is  summed  from  1 to  3 and  x3  = y.  c(y)  and  p(y) 
are,  assumed  to  be  Lebesgue  measurable  on  R+  = (y:  y > 0}  and  to 
satisfy 


0 < c:  < c(y)  < c2  < 


> for  all  y 6R+ 


0 < Pj  < p(y)  < p2  < ® i 

where  Cj,  c2,  pt  and  p2  are  suitable  constants. 


(7.7) 


7.3  Hilbert  space  formulation 
The  operator 


Au  = -c2(y)p(y) 


1 9u 


3x^  |p(y)  3xj 


(7.8) 
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was  shown  in  section  2 to  be  formally  selfadjoint  with  respect  to 
the  inner  product 


(u,v) 


u(x,y)  v(x,y) 


c-2 (y) 


p-1 (y)dxdy 


(7.9) 


where  dxdy  denotes  integration  with  respect  to  Lebesgue  measure 
on  R2.  The  corresponding  Hilbert  space  is 


= L2(R|,c-2(y)  p-1 (y)dxdy) 


(7.10) 


The  solution  of  the  initial-boundary  value  problem  (7.4)  - (7.6) 
given  below  is  based  on  a selfadjoint  realization  of  A in  3C.  To 
define  it  let  P(R+)  denote  the  Schwartz  space  of  R^.  and  P’ (R+) 
the  dual  space  of  all  distributions  on  R^..  The  Lebesgue  space 
L2(Rl)  can  be  regarded  as  a linear  subspace  of  P'(R+).  Note  that 
L2(R+)  and  If  are  equivalent  Hilbert  spaces  by  (7.7).  Let 


L*(R|)  = L2(Rj)  n {u:  au/Oxj  € L2(Rj),  j = 1,2,3}  (7.11) 


denote  the  first  Sobolev  space  of  R+.  It  is  a Hilbert  space  with 
inner  product 

3 

(u,v)j  = (u,v)0  + l Ou/3Xj ,3v/SXj)0  (7.12) 

J-l 

where  (u,v)0  is  the  inner  product  in  L2(R^).  P(R^_)  defines  a 

linear  subset  of  L2(R^)  and  hence 


L^’°(R*)  = closure  of  P(R2)  in  L*(R2) 


(7.13) 


is  a closed  linear  subspa  e of  L2(R^).  It  is  known  that  all  the 
functions  in  L2’°(R+)  satisfy  the  Dirichlet  boundary  condition 
(7.5)  as  elements  of  L2(R2);  see  [19]  and  [43, Cor.  2.7]. 

A realization  of  the  operator  A in  3f  will  be  defined  by 


D(A) 


Lj'°(R*)  n (u:  Au  e 3f> 


and 


(7.14) 


Au  *=  Au  for  all  u € D(A) 


(7.15) 
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To  interpret  the  condition  Au  € Jf  in  (7.14)  note  that  if  u 6 L*(j\p 
then  p-1(y)3u/3xj  € L2(r|)  for  j = 1,2,3.  The  second  derivative 
in  (7.8)  may  therefore  be  Interpreted  in  the  sense  of  V'  (R+)  . Thus 
Au  G V (R^)  and  the  condition  Au  G JC  is  meaningful. 

The  selfadjointness  of  A in  JC  may  be  proved  by  the  method  of 
[43, §2].  Another  proof  may  be  based  on  the  theory  of  sesquilinear 
forms  in  Hilbert  space  [18, Ch.  6].  These  methods  imply  the 


7.4  Theorem 


A is  a selfadjoir.t  real  positive  operator  in  JC.  Every  u G D(A) 
satisfies  the  Dirichlet  boundary  condition  (7.5)  as  an  element  of 
L2(R2).  Moreover,  D(A^2)  = L2>°(R^)  and 

3 

BA^ull2  = l 
j=l 


| 3u/9xj | 2P-1 (y)dxdy 


for  all  u e D(A1/2) 


(7.16) 


The  operator  A may  be  used  to  construct  "solutions  in  JC'  and 
"solutions  with  finite  energy"  of  (7.4)  - (7.6),  as  described  in 
section  2.  The  detailed  analysis  of  the  structure  of  these  solu- 
tions will  again  depend  on  the  construction  of  an  eigenfunction 
expansion  for  A.  For  simplicity,  the  construction  will  be  described 
here  for  a special  choice  of  the  functions  c(y)  and  p(y).  Never- 
theless, the  results  obtained  are  typical  of  a large  class  of 
stratified  fluids. 


7.5  The  Pekeris  model 

This  name  will  be  used  for  the  stratified  fluid  defined  by 


c(y)  = < 


P(y)  = < 


Cl, 

0 

< 

y < h 

C2  » 

y 

> 

h 

Pi* 

0 

< 

y < h 

P2  » 

y 

> 

h 

1 » 

p2  and 

: h are 

(7.17) 


(7.18) 


sitive  constants.  This  model  was 
used  by  C.  L.  Pekeris  [26]  in  his  study  of  acoustic  wave  propaga- 
tion in  shallow  water.  The  model  represents  a layer  of  water  with 
depth  h,  sound  speed  and  density  pt  which  overlays  a bottom, 

such  as  sand  or  mud,  with  sound  speed  c2  and  density  p2 . 


Water 


X - (X!,X2) 


Sound  Speed  Cj 
Density  pj 


Wet  Sand 


Sound  Speed  c2 
Density  p2 


Figure  7.  The  Pekeris  Model 


The  most  interesting  case  occurs  when 


Cj  < c2 


(7.19) 


and  this  condition  is  assumed  to  be  satisfied  in  what  follows.  A 
detailed  study  of  the  Pekeris  operator  was  given  by  the  author  in 
[43] . Here  the  main  results  of  [43]  are  reviewed  and  used  to 
derive  the  asymptotic  wave  functions  for  the  Pekeris  model. 


7.6  Eigenfunctions  of  A 

It  was  shown  in  [43]  that  A has  a pure  continuous  spectrum  r.nd 
a complete  family  of  generalized  eigenfunctions  was  constructed. 
These  functions  w(x,y)  are  characterized  by  the  following 
properties 


w is  locally  in  D(A) 

Aw  = Aw  for  some  A > 0 

w(x,y)  is  bounded  in  R+ 

w(x.v)  « (2tt)_1  eip  x w(y),  p € R2 


(7.20) 

(7.21) 

(7.22) 
(7.2V, 


where,  in  (7.23),  w(y)  is  independent  of  x.  The  eigenfunctions  are 
of  two  types,  called  free  wave  eigenfunctions  and  guided  wave 
eigenfunctions.  Their  definitions  and  physical  interpretations 
follow. 
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7.7  Free  wave  eigenfunctions 

These  functions  exist  when  the  eigenvalue  A satisfies 


A > c| Jp j 2 > c2|P|2,  j p j 2 = p^  + p^ 


(7.24) 


To  define  them  let 


C = (A/c|  - | p | 2) 1/2  >0,  n = (A/c2  - |P | 2) 1/2  > 0 (7.25) 


w0(x,y,p,A)  = (2tt)-1  elp  x w (y,p,A) 


(7.26) 


where 


sin  riy 


wo (y>P>*)  “ a(p,A) < 


, 0 < y < h 


(7.27) 


Y+(C.n)  eiC(y_h)+Y_(4,n)e~iC(y"h),  y > h 


Y+(C»n)  = \ sill  nh  + ^ ^ cos  nh 


(7.28) 


In  (7.27),  a(p,A)  is  a positive  normalizing  constant.  It  was 
shown  in  [43]  that  the  eigenfunction  expansion  takes  its  simplest 
form  when 


a(p,A)  = p^2/2(ttO  1/2  |y+(C,H) 


(7.29) 


In  physical  terms,  the  eigenfunction  w0(x,y,p,A)  represents 
an  acoustic  field  with  time  dependence  exp  (-itA1'2)  which  is  the 
sum  of  two  plane  waves  in  each  layer.  It  may  be  interpreted  as  a 
plane  wave  which  propagates  in  the  region  y > h,  is  refracted  at 
y = h,  reflected  at  y = 0 and  refracted  again  at  y = h;  see 
Figure  8 where  the  propagation  directions  are  indicated.  It  can 
be  verified  that  Snell's  law  of  refraction  is  satisfied  at  y = h 
and  the  law  of  reflection  is  satisfied  at  y = 0. 
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Figure  8.  Ray  diagram  for  free  wave  eigenfunction 


7.8  The  dispersion  relation 

For  values  of  X which  satisfy 

ci  I P 1 2 < A c | p | 2 (7.30) 

the  function  w0(x,y,p,A)  defined  by  (7.25)  - (7.28)  still  satis- 
fies conditions  (7.20),  (7.21)  and  (7.23).  However,  (7.30)  im- 
plies that  £ is  pure  imaginary,  say 

i - V - ( |p  | 2 - A/c*)l/2  > 0 (7.31) 

while  r|  is  still  real  and  positive.  It  follows  that  w0(x,y,p,A) 
satisfies  the  boundedness  condition  (7.22)  if  and  only  if 

Y_(i£' ,n)  = 0 (7.32) 

or,  by  (7.28), 


V 


_ £2. 
Pi 


ri  ctn  nh 


(7.33) 


For  X and  |p|  which  satisfy  (7.30),  (7.33)  is  equivalent  to  the 
sequence  of  equations 
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hn  = ( k — ~)  7r  + tan-1 


PiC 


P2n 


, k = 1,2, 


(7.34) 


where  |tan-1  a|  < tt/2.  Each  equation  (7.34)  defines  a functional 
relation  between  |p|  and  X or,  equivalently,  between  |p|  and 


a)  = Xl/z 

The  solutions,  which  will  be  denoted  by 

X = \(|p|).  = u»k(  |p|)  = xk(  I P I ) 1/2 


(7.35) 


(7.36) 


represent  a relation  between  the  wave  number  |p]  of  the  plane 
waves  in  w0(x,y,p,X)  and  the  corresponding  frequencies  U).  Such 
relations  are  called  dispersion  relations  in  the  theory  of  wave 
motion.  The  relations  (7.34),  (7.36)  were  analyzed  in  [43]  and 
found  to  have  the  following  properties. 


7.9  Properties  of  co^ ( | p | ) 

For  each  k = 1,2,3,  •••  define 


Then 


Pi  =»  TtCj /2h(c2  - c2)l/z, 

pk  = <2k  " 1)pi 

(7.37) 

w^(|p|)  is  analytic  and 

^k ( 1 ? I > > 0 for  IpI  _ Pk 

(7.38) 

ci Ip!  5 “k< 1 P 1 > < c2 1 P 1 

|p|  > Pk 

(7.39) 

wk(Pk>  = c2Pk’  = 

C2 

(7.40) 

a)k(  1 p 1 > ~ ci  IpI  for  lp  1 

->  00 

(7.41) 

Moreover,  an  explicit  parametric  representation  of  the  dispersion 
curves  (7.36)  was  given  in  [43]. 


7.10  Guided  wave  eigenfunctions 

The  functions  w^x.y.p)  = w0  (x,y, p,X^( | p | ) ) are,  by  construc- 
tion, the  solutions  of  (7.20)  - (7.23)  for  eigenvalues  which 
satisfy  (7.30).  It  was  shown  in  [43]  that  there  are  no  solutions 
of  (7.20)  - (7.23)  when  X < c2|p|2.  The  functions  w^(x,y,p)  have 
the  form 
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wk(x,y,p)  = (2tt)~  1 eip  x wk(y,p) 


(7.42) 


where 

sin  nk(|pj)y  , 0 < y < h 

wk(y.P>  = a,  (p)  < (7.43) 

, , -q(|P|)(y-h) 

sin  nk(|p|)h  e , y > h 

with 

\(|p|)  = (\(|p!)/Cj  - I p | 2) 1/2 , 

> (7.44) 

Ck(|p|)  = ( I P 1 2 - Ak(|p|)/c|)1/2 

In  (7.43),  ak(p)  is  a positive  constant  which  is  determined  by 
the  condition 
>00 

|wk(y*P)|2  c-2(y)  p-1 (y)dy  = 1 (7.45) 

« o 

In  physical  terms,  the  eigenfunction  wk(x,y,p)  represents  an 
acoustic  field  with  time  dependence  exp  (— i tcok ( | p | ) ) which  corre- 
sponds to  a plane  wave  which  is  trapped  in  the  layer  0 < y < h by 
reflection  at  y = 0 and  total  internal  reflection  at  the  interface 
y = h.  In  the  layer  y > h the  field  is  exponentially  damped  in 
the  y-direction  and  propagates  strictly  in  the  horizontal  direction 
p;  see  Figure  9 where  the  propagation  directions  are  indicated. 


Air 


Exponential  — > — — > 

Damping 


> — — * 

Figure  9.  Ray  diagram  for  guided  wave  eigenfunctions 
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7.11  The  eigenfunction  expansion 

The  free  wave  eigenfunctions  w0(x,y,p,A)  are  parameterized  by 
the  region 


G0  = {(p,A):  p e R2  and  c2|p|2  < A}  C R3 


(7.46) 


Similarly,  the  guided  wave  eigenfunctions  w^(x,y,p)  are  parameter- 
ized by  the  regions 


= tP:  I P I > Pkh  k = 1,2,' 


(7.47) 


The  following  expansion  theorem  was  proved  in  [43],  First  of  all, 
for  every  f <=  the  limits 


Mp.*) 


L2(G0)-lim 

M-*c°  J 


M 


(7.48) 


w0(x,y,p,A)  f(x,y)c"2(y)p  l(y)dxdy 
I x I <M 


and 


Mp)  = 


L2  (Gk)-lim 
M-x» 


M 


(7.49) 


wk(x,y,p)  f(x,y)c“2(y)p  '(yjdxdy, 
I x I <M 


k = 1,2, 


exist  and  satisfy  the  Parseval  relation 


,f“^  = lfk"L2(nk) 


(7.50) 


Moreover,  if 


G^  *=  { (p, A) : p e R2  and  c|  { p l 2 < A < M> 
G^  = ^P:  Pk  < I P I <M>,  k = 1 , 2, • • • 


(7.51) 


then  the  limits 
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f0(x»y)  “ JC-lim  w0(x,y,p,A)  f0(p,A)dpdA  (7.52) 

M-x«  J„M 


fk(x,y)  = 3C-lim  f w (x,y,p)  f,(p)dp,  k = 1,2, •••  (7.53) 

k 

exists  and  satisfy 


f(x,y)  = JC-1  im  £ f,(x,y) 


(7.50 


The  relations  (7. 48),  (7.49),  (7.52),  (7.53)  and  (7.54)  will  also 
be  written  in  the  following  more  concise  symbolic  forms,  in  analo- 
gy with  the  notation  of  previous  sections. 


f0(P>X)  = w (x,y,p,A)  f (x,y)  c-2(y)  p-1(y)dxdy  (7.55) 


fk(p)  = wk(x,y,p)  f (x,y)  c-2 (y)  p'^yjdxdy, 


k = 1,2,' 


Mx*y>  “ w0 (x , y , p , A ) f 0 (p, A)dpdA 


fk(x,y) 


wk(x,y,p)  fk(p)dp,  k = 1,2,*' 


(7.56) 


(7.57) 


(7.58) 


f(x,y)  = l fk(x,y) 


(7.59) 


Equations  (7.55)  - (7.59)  are  the  eigenfunction  expansion  for 
A and  show  the  completeness  of  the  generalized  eigenfunctions  de- 
fined above.  The  representation  is  a spectral  representation  for 
A in  the  sense  that,  for  every  f € D(A) , 

(Af)0(p,A)  = A f0(p,A) 


(7.60) 
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and 


(Af)k(p)  = Ak(|p|)  fk(p),  k = 1,2, ••• 


(7.61) 


The  representation  (7.55)  - (7.59)  defines  a modal  decomposi- 
tion for  the  Pekeris  model.  It  was  shown  in  [43]  that  if 

\ = (fk:  f € X]  c X,  k = 0,1,2,  •••  (7.62) 

then  each  X is  a closed  subspace,  3C  and  X.  are  orthogonal  for 
k t Z and  k x. 


JC=  I © \ (7.63) 

k=0 


Moreover,  it  was  shown  that  (7.60),  (7.62)  imply  that  (7.63)  re- 
duces A.  In  fact,  more  was  shown  in  [43];  namely  that 

$Rf  = fk  e L2(£y,  k = 0,1,2,  •••  (7.64) 

defines  an  operator 

X + L2(flk),  k = 0,1,2,  •••  (7.65) 

which  is  a partial  isometry  with  initial  set  JC  and  final  set 
L2(ftk);  i.e., 

K \ = V \ $k  * 1> k = t*’1’2’"'  (7.66) 

where  is  the  orthogonal  projection  of  X onto  3^. 


7.12  Solution  in  3f  of  the  propagation  problem 

Attention  will  again  be  restricted  to  the  case  where  f £ 3f 
and  g e D(A-1^2)  so  that  the  solution  in  3C  has  the  form 

u(t,x,y)  = Re  (v(t,x,y)}  (7.67) 


with 


v(t,  • , •)  = exp  (-itA1,/?)h,  h = f + iA~1,/2g  e 3f 


(7.68) 
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The  modal  decomposition  of  v(t,x,y)  is 


'(t,x,y)  = l v,  (t,x,y) 


k=0 


where 


(7.69) 


v0(t,x,y)  = 


w (x,y,p, A)  exp  (-itA1/2)  h (p,A)dpdA  (7.70) 


and 


vk(t,x,y)  = 


wk(x,y,p)  exp  (-it(ak(|p|))hk(p)dp, 


n,. 


k = 1,2,' 


(7.71) 


Moreover,  the  modal  waves  vk(t,x,y)  are  independent  in  the  sense 
that  they  are  orthogonal  in  JC  for  every  t 6 R because  (7.63)  is  a 
reducing  decomposition  of  A.  Asymptotic  wave  functions  for  each 
mode  will  now  be  calculated  beginning  with  the  guided  modes  v,  , 
k > 1.  * 


7.13  Asymptotic  wave  functions  for  the  guided  modes  (k  > 1) 


If  the  representation  (7.42)  for  the  eigenfunctions  wk(x,y,p) 
is  substituted  into  (7.71)  the  spectral  integrals  takes  the  form 

vk(t,x,y)  = 

(7.72) 

exp  {i(x*p-twk(  |p|))}wk(y,p)hk(p)dp,  k = 1,2,“* 


JL_ 

2tt 


where  wk(y,p)  is  defined  by  (7.43).  The  behavior  for  large  t of 
these  integrals  will  be  calculated  by  the  method  of  stationary 
phase.  In  the  present  case  the  integral  is  a double  integral 
(^k  C R2)  and  the  phase  function 

0k(p,x,t)  = x*  p - twk(|p|)  (7.73) 


is  stationary  with  respect  to  p if  and  only  if 


In  particular,  the  number  and  distribution  of  the  stationary 
points  is  determined  by  the  group  speed  function  for  the  kth 
guided  mode: 

\(|p|)  = “£(|p|)»  I P I > Pk  (7.7 

The  defining  relation  (7.34)  for  w^(|p|)  implies  the  following 

7.14  Properties  of  U^(|p|) 

For  each  k = 1,2,3,* ••  there  exists  a unique  p,  > p,  where 
uk(pk)  = 0 and  pR  > pR.  Moreover, 

0 < Uk  = Uk(p^)  < Uk(|p|)  < c2  for  all  |p|  > pR  (7.7 

Uk(|p|)  < 0 for  pk  < |p|  < pk  and  U^(|p|)  >0 

for  J p | > pk 

lim  uk( I P I ) = c2>  lim  uk(|pl)  = ci  (7-7i 

|p|-*pk  Ip  I-*00 


These  properties  are  indicated  in  Figure  10. 


The  stationary  points  of  0k(p,x,t)  are  defined  by  (7.74). 
This  may  be  written  in  (2-diinensional)  vector  notation 
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x/t  = Uk(|p|)p/|p|  (7.79) 

This  is  equivalent  to  the  conditions 

Uk(|p|)  = |x|/t  (7.80) 

and 


p is  parallel  to  x and  in  the  same  direction 


(7.81) 


s kice  Uk ( | p | ) > 0 and  t is  assumed  to  be  positive.  Conditions 
(7.80)  and  (7.81)  determine  |p|  and  p/|p|,  respectively.  In 
particular,  it  is  clear  from  Figure  10  that 


For  | x I > c2t  and 
stationary  phase 


i i A 

|x|  < Ukt  there  are  no  points  of 


For  Cit  < |x|  < c2t  and  |x|  = U^t  there  is  one  point  g2^ 
of  stationary  phase 


For  uft  < | x I < Cjt  there  are  two  points  of 
stationary  phase 


According  to  the  method  of  stationary  phase  each  stationary 
point  where  det  (820k/3pi3p.)  ^ 0 (regular  stationary  point)  con- 
tributes a term 

exp  (i0,  (p,x,t)  + i 7 sgn  (320, /3p.3p .)) 

£ 2 ^ Wk(y,P)hk(p)  (7.83) 

| det  (320k/3Pi3Pj) | 1/2 

to  the  asymptotic  expansion  of  the  integral  (7.72),  where  sgn  and 
det  denote  the  signature  and  determinant,  respectively,  of  the 
Gramian  matrix  (320k/3p .3p .) . A short  calculation  shows  that  the 
eigenvalues  of  the  GramianJfor  (7.73)  are  -tUk(|p|)  and  -tUk(|p|)/ 
and  hence  for  t > 0 

sgn  (320k/3Pl3pj)  = -1  - sgn  Uk(|p|)  (7.84) 

det  (320k/3p.3p.)  = t2Uk( | p| )Uk( 1 P | ) / | p | 


(7.85) 
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In  particular,  the  stationary  points  are  regular  when  |p|  jf  p.. 
Substitution  of  (7.73),  (7.84)  and  (7.85)  into  (7.83)  gives 
the  function 


vt(t,x,y,p)  = 


|p|l/2  exp  (i(|x|  ] p | - tcok ( | p | ) -f  - J sgnU^(|p|))} 

t f Uk ( I P I ) |U’ ( |p|)> 1/2 


(7.86) 

wk(y.p)\(p) 


To  find  the  asymptotic  wave  function  for  vk(t,x,y)  it  is  necessary 
to  solve  (7.79)  for  p and  substitute  in  (7.86).  The  result  may  be 
described  by  means  of  the  two  inverse  functions  of  U^(|p|)  which 
may  be  defined  as  follows:  see  Figure  10. 


Ip  I “ 0 uk(|p|)  = q and  p.  < |p|  < p£ 


! P I " pk(q)  ° U.(|p|)  = q and  |p|  > p£ 


(7.87) 


It  is  clear  from  the  discussion  of  Uk(|p|)  that  Pk  and  P[<  are  analy- 
tic functions,  Pk  maps  {q:  < q < c/}  onto  { j p| : pk  < |p)  < p£} 

and  Pk  maps  { q : < q < c,}  onto  Tip]:  |p|  > p^} . 

The  asymptotic  behavior  of  vk(t,x,y)  can  now  be  described. 

The  point  of  stationary  phase  |p|  = p£(jx|/t)  makes  a contribution 

vr,f  (:.,x,y)  = Xu(t,x)  v!”(t  ,x,y  ,p£  ( |x  | /t)x/ |x  | ) (7.88) 


where 


X^(t,x)  is  the  characteristic  function  of 


{ ( t , x ) : U,  < | x | / 1 < c2) 


(7.89) 


Similarly,  the  point  of  stationary  phase  |p|  = P^(|x|/t)  makes  a 
contribution 


v”’S(t,x,y)  = xf(t.x)  v"(t,x,y,P*( |x|/t)x/ jx|) 


(7.90) 


where 
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Xjc(t,x)  is  the  characteristic  function  of 

(7.91) 

{(t,x):  U£  < |x|/t  < Cj} 

00  f 

The  functions  v^’  are  called  the  "fast  waves"  because  they  de- 
scribe waves  which  arrive  at  points  (x,y)  at  times  t = |x|/c2 
corresponding  to  the  speed  c(y)  = c2  of  waves  in  the  "fast"  me- 
dium filling  y > h.  Similarly,  the  functions  v^’s  are  called  the 
"slow  waves"  because  they  describe  waves  which  arrive  at  (x,y)  at 
times  t = |x|/cj  corresponding  to  the  speed  c(y)  = Cj  of  waves  in 
the  "slow"  medium  filling  0 < y < h.  Finally,  the  total  asymptotic 
wave  function  is  the  sum 

vk(t,x»y)  = v^’f(t,x,y)  + v“,s(t,x,y)  (7.92) 

The  following  convergence  theorem  was  proved  in  [40]  by  the  method 
outlined  in  section  5. 


7.15  Theorem 

Let  h e JC.  Then  for  each  k > 1,  vjc(t,*,*)  € 3f  for  all  t > 0 
and  t -*■  vjc(t,*,*)  € 3C  is  continuous.  Moreover,  v^(t,*,*)  is  an 
asymptotic  wave  function  for  the  modal  wave  vjc(t,*,*) 

■ exp  (-itA1^2)  P^h,  i.e., 


lira  I vR  ( t , • , • ) - v.  (t,*,  Olj, 

t*+oo 


= 0 


(7.93) 


The  same  methods  were  used  in  [40]  to  prove  convergence  in 
the  energy  norm  when  h has  finite  energy. 


Note  that  v^(t,x,y)  represents  a guided  wave  which  propagates 
radially  outward  in  horizontal  planes  y = const,  and  is  exponential- 
ly damped  in  the  vertical  coordinate  y.  This  is  evident  from  the 
defining  equations  (7.86),  (7.88),  (7.90)  and  (7.92). 


7.16  Asymptotic  wave  functions  for  the  free  mode 


It  will  now  be  shown  that  the  free  mode  wave  function  vQ(t,x,y) 
is  asymptotically  equal  in  If  to  a free  wave  propagating  with  speed 
c2  in  the  half-space  y > h.  To  this  end  note  that 


% 


{ i(x* p - tAl/2)}  w0 (y,p,A)h0 (p,X)dpdX 


v0(t,x,y)  = ~ 


exp 


(7.94) 
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where  h0  € L2(ft0).  The  representation  (7.27)  for  w0(y,p,A) 
implies  that 


v0(t,x,y)  = v (t,x,y-h)  + v“(t,x,y-h),  y > h 


(7.95) 


where 


v0(t,x,y)  = 


i(x*p+y£;-tA1/2) 


a(p,A)  Y+(5»n)h0(p,A)dpdA 
(7.96) 


v0(t,x,y)  = ■—  ex^x  P y?  tA  * a(p,A)  y_(£;,n)h0  (p,  A)dpdA 


(7.97) 


The  change  of  variables 


(p,A)  **•  (p»q)»  q = C = (A/c2  - | p | 2 ) 1/2 


(7.98) 


in  (7.96)  gives 


v„ (t,x,y)  = 


(27T)  3/2 


i(x*p+yq-tOj(p,q/) 


h(p»q) dpdq 


(7.99) 


where 


h(p,q)  = c2p2l/2(2|q|)1/2(y+(5,n)/jY+(C,n)  |)h0(p,A)  (7.100) 


A = A(p,q)  = to(p,q)2  = c2(|p|2  + q2) 


Similarly,  the  change  of  variables 


(p,*)  (p.q).  q = = -(X/c|  - | p | 2 ) 1/2 


(7.101) 


(7.102) 


in  (7.97)  gives 


i (x  • p+y  q-  to)  (p , q)  ) £ 


h(p,q)dpdq 


(7.103) 


where  h is  defined  by  (7.100).  Adding  (7.99)  and  (7.103)  and 
using  (7.95)  shows  that 
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v„(t,x,y+h) 


1 

(2rr)  3/2 


ei(x.p+yq-ta)(p,q))fi(p>q)dpdq 

(7.104) 


for  all  y > 0.  Moreover,  (7.100)  implies  that  h G L2(R3) • Thus 
(7.104)  and  (7.101)  imply  that  in  the  half-space  y > h vn(t,x,y) 
coincides  with  a solution  in  L2(R3)  of  the  d'Alembert  equation 
with  propagation  speed  c2.  Now,  the  results  of  section  3 imply 
that  the  right-hand  side  of  (7.104)  has  an  asymptotic  wave  function 
in  L2 (R3) ; say 


w (t,x,y)  = G(r- c2t,0)/r,  r2  = |x|2  + y2,  G = (x,y)/r 


It  follows  that  if 


(7.105) 


oo  f w (t,x,y-h)  , y > h 

v0(t,x,y)  = < 

^ 0,  0 < y < h 


(7.106) 


then 


lim  |v0(t,*,  •)  - v^t,  • , *)ljf  = 0 (7.107) 

£—►00 

A proof  may  be  found  in  [40]  . This  paper  also  contains  a proof  of 
convergence  in  the  energy  norm,  when  h has  finite  energy,  and 
applications  of  these  results  to  the  calculation  of  asymptotic 
energy  distributions  in  stratified  fluids. 


7.17  Other  cases 

The  case  of  the  symmetric  Epstein  profile,  defined  by 

c“2(y)  = c~2  sech2(y/H)  + c~2tanh2 (y/H)  (7.108) 

and  p(y)  = 1 was  studied  by  the  author  in  [41]  where  eigenfunction 
expansions  and  asymptotic  wave  functions  are  derived.  Eigenfunc- 
tion expansions  for  the  case  of  the  general  Epstein  profile 


c~2(y)  = K cosh2 (y/H)  + L tanh  (y/H)  + M 


(7.109) 
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and  p(y)  = 1 have  been  given  by  J.  C.  Guillot  and  the  author  [1A, 
15] . Asymptotic  wave  functions  fcr  this  case  are  currently  being 
studied  by  Y.  Dermenjian,  J.  C.  Guillot  and  the  author.  Prelim- 
inary studies  show  that  the  results  given  above  for  the  Pekeris 
model  are  valid  for  a large  class  of  profiles  c(y),  p(y) . The 
essential  hypotheses,  apart  from  the  boundedness  (7.7),  are  that 
c(y)  should  have  a global  minimum  at  some  finite  point  and  that 
c(y)  should  tend  to  a limit  at  infinity  sufficiently  rapidly,  if 
c(y)  does  not  have  a minimum  then  there  are  no  guided  waves. 
However,  these  results  have  not  yet  been  proved  in  this  generality. 


8.  PROPAGATION  IN  CRYSTALS 

Acoustic  wave  propagation  in  an  unlimited  homogeneous  crystal 
is  analyzed  in  this  section.  The  analysis  is  similar  to  that  for 
homogeneous  fluids  given  in  section  3.  The  principal  new  feature 
is  the  influence  of  anisotropy  on  the  structure  of  the  asymptotic 
wave  functions. 

A homogeneous  crystal  is  characterized  by  a constant  density 
p(x)  = p and  stress-strain  tensor  cjj^(x)  = c-jjjj.  It  will  suffice 
to  consider  the  case  p = 1.  Thus  the  propagation  problem  reduces 
in  this  case  to  the  Cauchy  problem  for  the  system 


9S  £m  9S 

J-  _ P 


3t‘ 


'jk 


9xk9xn 


= 0,  j = 1,2,3 


where  the  constants  c^f1  satisfy  (2.13)  and  (2.35). 

J k 


(8.1) 


8.1  Hilbert  space  formulation 

It  was  shown  in  section  2 that  the  differential  operator  A 
defined  by 


(A„,j 


£m 

Cjk 


a2un 


9x.  9x 
k m 


j 


1,2,3 


is  forma] ly  selfadjoint  in  the  Hilbert  space  K 
inner  product 


(8.2) 


L2(R3.C’)  with 


1 R 


3 


(u,v) 


Uj  (x)  Vj  (x)  dx 


(8.3) 
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In  fact,  the  operator  A in  X with  domain  D(A)  = 0(R3)  is  essen- 
tially selfadjoint  and  its  unique  selfadjoint  exten  :i  is  the 
operator  A defined  by 

D(A)  = X n {u:  Au  e 3f}  (8.4) 

Au  = Au  for  all  u £ D(A)  (8.5) 

It  is  easy  to  verify,  using  the  Plancherel  theory  of  the  Fourier 
transform,  the  following 

8.2  Theorem 

A is  a selfadjoint,  real  positive  operator  in  X. 

It  follows,  as  in  preceding  sections,  that  the  Cauchy  problem 
for  (8.1)  has  a solution  in  X of  the  form 

u j ( t , x ) = Re  (vj (t ,x) } (8.6) 

where 

v(t,‘)  = exp  (-itA^2)h,  h = f + iA-1|/2g£  X (8.7) 

whenever  the  Cauchy  data  u(0,x)  = f(x)  and  3u(0,x)/9t  = g(x)  satisfy 
f £ X,  g £ D(A_1/2) . 

8.3  Fourier  analysis  of  A 

The  Plancherel  theory  of  the  Fourier  transform  $0:  L2 (R3) ** L2 (R3) 
was  defined  and  used  in  section  3;  see  (3.12).  It  may  be  extended 
immediately  to  X = L2(R3,C3)  by  defining 

$>0U  = *0(ui,u2,u3)  = ($0u1  ,$0u2,$0u3)  (8.8) 

and  <i>0  is  also  unitary  in  X.  Property  (3.14)  implies  that  the 
operator  $0A  4>*  corresponds  to  multiplication  by  the  3 * 3 matrix 
valued  function 

A(p)  = (Aj£(p))  = (Cj“  PkPm).  P e R3  (8-9) 


Thus 
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A = <*>*  A( • ) <J>Q  (8.10) 

Moreover,  conditions  (2.13)  and  (2.35)  Imply  that  A(p)  is  a real 
Hermitian  positive  definite  matrix  for  all  p €.  R3  - {0}.  The 
spectral  analysis  of  A will  be  based  on  (8.10).  The  analysis 
begins  with 


, 


8.4  Spectral  analysis  of  A(p) 

The  eigenvalues  of  A(p)  are  the  roots  p of  the  characteristic 
polynomial 

det  (pi  - A(p) ) = 0 (8.11) 


The  Hermitian  positive  definiteness  of  A(p)  implies  that  the  roots 
are  real  and  positive  for  all  p E R3  - {0} . They  may  be  uniquely 
defined  as  functions  of  p ’ey  enumerating  them  according  to  their 
magnitudes: 

0 < Pj (p)  < p2(p)  < y3(p)  for  all  p e R3  (8.12) 

A result  of  T.  Kato  [18]  implies  (see  also  [37]) 

p j : R3  ■*  R is  continuous,  j = 1,2,3  (8.13) 

Equation  (8.11)  implies  that  p^ (p)  is  homogeneous  of  degree  2 

Pj(ap)  = ct2Pj(p)  for  all  a 6 R and  p £ R3  (8.14) 

The  functions 


A^p)  = /pj(p)',  p € R3,  j = 1,2,3 


(8.15) 


are  also  needed  below.  A detailed  study  of  these  functions  has 
been  made  by  the  author  in  connection  with  a formulation  of  elas- 
ticity theory  in  terms  of  first  order  symmetric  hyperbolic  systems; 
see  [29,35,37,44].  A number  of  results  from  these  papers  are 
quoted  and  used  below. 


It  was 
nomial  0(p) 
roots  p . (p) 


shown  in  [44]  that  there  exists  a homogeneous  poly- 
1 0 such  that  the  points  p e R3  where  two  or  more 
coincide  are  contained  in  the  cone 


ii 


I 


J 
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Z = {p  e R3:  0(p)  = 0} 


(8.16) 


Thus 


o < yj(p)  < U2(p)  < P3(p)  for  all  p € R3  - Z (8.17) 

It  follows  that 

Pj (p)  is  analytic  on  R3  - Z,  j = 1,2,3  (8.18) 

The  orthogonal  projection  of  C3  onto  the  eigenspace  for  y.(p)  is 
given  by  [18]  3 


(A(p)  - zl)_1dz,  j = 1,2,3  (8.19) 

VP) 

where 

Yj (p)  = {z:  |z  - Mj (p) | = Cj(p)},  j = 1,2,3  (8.20) 

and  the  radii  c^(p)  are  chosen  so  small  that  the  3 circles  Yj(p) 
are  disjoint.  This  is  possible  for  all  p £ R3  - Z by  (8.17).  The 
matrix  valued  functions  so  defined  can  be  shown  to  have  the 
following  properties  [18,44]: 


Pj(p)  is  analytic  on  R3  - Z,  j = 1,2,3  (8.21) 

Pj(ap)  = Pj(p)  for  all  a ft  0 (8.22) 

Pj(p>*  = VP)’  VP)  Pk(p)  = 6jkPk(p)  for  P6r3"z  <8-23) 

3 

l Pj(p)  = 1 for  p £ R3  - Z (8.24) 

j=l 

A(p)  Pj(p)  = M j (p)  Pj(p)  for  p € R3  - Z,  j = 1,2,3  (8.25) 


The  last  two  properties  imply  that  the  projections  Pj(p)  de- 
fine a spectral  representation  for  A(p) ; i.e., 


J 
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3 

A(p)  = l ^(p)  Pj(p)  for  p e R3  - Z (8.26) 

j=l 

8.5  Spectral  analysis  of  A 

The  representations  (8.10)  and  (8.26)  provide  a complete 
spectral  analysis  of  A.  In  particular,  it  follows  that  A is  an 
absolutely  continuous  operator  whose  spectrum  is  [0,<»)  (cf.  [36, 
44]).  Moreover,  if  ^((j)  is  any  bounded  Lebesgue-measurable  func- 
tion of  ii  > 0 then 

3 

'J'(A)  = 9*  l f(bk(*))  Pfc(*)  4>0  (8.27) 

k=l 


8.6  Solution  in  3f  of  the  Cauchy  problem 

Application  of  (8.27)  to  the  solution  in  3C  (8.7)  yields  the 
representation 


v(t,x)  = l vk(t,x) 


(8.28) 


k-1 


where 


vk(t’x)  = 


i(x*p-tXk(p)) 


J R3 


Pk(p)h(p)dp  (8.29) 


and  Xk(p)  = /uk(p) . Of  course,  the  integral  in  (8.29)  converges 
in  3f,  in  the  sense  of  the  Plancherel  theory,  rather  than  pointwise. 
Equations  (8.28),  (8.29)  represent  solutions  in  3f  of  (8.1)  as  a 
superposition  of  solutions 

i(x*p-tX. (p))  ^ 

e K Pk(p)h(p)  (8.30) 


This  may  be  interpreted  as  a plane  wave  which  propagates  in  the 
crystal  with  direction  p/|p|,  wave  number  |p|  and  frequency 


to  = Xk(p) 


(8.31) 


81 


The  polarization  of  the  wave  is  determined  by  P^(p) • The  corres- 
ponding generalized  eigenfunctions  of  Z'  are  the  matrix  plane  waves 
[30] 


wk(x,p) 


1 

(2 IT)  3/2 


exp  (ix*p)  PR(p) 


(8.32) 


8.7  The  dispersion  relation,  phase  and  group  velocities 

The  dispersion  relation  between  the  frequency  oj  and  wave  vec- 
tor p of  plane  waves  in  the  crystal  is  (8.31)  or,  by  (8.15)  and 
(8.11) 

det  (o>2l  - A(p) ) = 0 (8.33) 


The  phase  velocity  for  (8.1)  is 

A,  (p)  , , 

vph('»  ■ M • i£r ' jP7T p ■ (itr)  ifr  (8 

by  the  homogeneity  of  A^(p) . The  group  velocity  for  (8.1)  is 


vg(p)  = Vp<°  = VpXk(p)  (8.35) 

The  medium  is  said  to  be  isotropic  if  v fl(p)  and  v (p)  have  the 
same  direction  for  all  p £ R3  - {0}.  Otherwise  it8is  said  to  be 
anisotropic.  It  is  easy  to  verify  that  the  medium  is  isotropic 
if  and  only  if  A (p)  is  a function  of  |pl  alone.  In  this  case 

\(p)  = ClJpl  ana  vph(p)  = vg(p)  = ckp/'pl* 

The  phase  and  group  speeds  for  (8.1)  are  the  magnitudes  of 
the  corresponding  velocities.  Thus 

Cph(p)  = lVph(P>l  = Ak(p/|pl)  1 

> (8.36) 

cg(p)  = lvg(p)l  = VpAk(P)  J 

Note  that  both  are  homogeneous  of  degree  zero  in  p and  hence  depend 
only  on  the  direction  of  propagation  p/|p|.  The  anisotropy  of  the 
medium  characterized  by  (8.1)  can  be  visualized  by  means  of 
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8.8  The  slowness  surface  S 

This  Is  the  real  algebraic  variety  defined  by 

S = {p  e R3:  det  (1  - A(p) ) = 0}  (8.37) 

It  is  clear  from  the  definition  of  the  Ak(p)  that 

■-  U 

k=l 

where 

Sk  = {p  e R3: 

or,  by  (8.36)  and  the  homogeneity  of  A,  (p) , 

Sk  = (P  e R3:  |p|ck(p)  = 1)  (8.40) 

Thus  p £ S if  and  only  if  |p|  is  the  reciprocal  of  a phase  speed 
for  the  direction  p.  Note  that  the  slowness  surface  of  an  isotropic 
medium  is  a set  of  concentric  spheres  with  centers  at  the  origin. 

The  properties  of  slowness  surfaces  were  studied  in  [35]  and 
[ 44 ] . In  particular,  the  following  properties  were  established 

Sk  is  continuous  and  star-shaped  with  respect  to  0 (8.41) 

As  an  algebraic  variety,  S will  in  general  have  singular  points 
and  these  are  precisely  the  set 

Zl  = (p  € S:  o £ S,  n S.  for  some  j ^ k)  (8.42) 

b ■ j k. 

Hence 

Sk  - Zg , k = 1,2,3,  are  disjoint  and  analytic  (8.43) 


8.9  The  wave  surface  W 

The  variation  of  the  phase  speed  with  direction  is  represented 
by  the  slowness  surface  S.  Similarly,  the  variation  of  the  group 
speed  is  represented  by  the  wave  surface  W.  W may  be  defined  as 


J 


Ak(p)  = 1} 


(8.38) 


(8.39) 
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the  polar  reciprocal  of  S with  respect  to  the  unit  sphere.  This 
means  that 

W = {x  6 R3:  x • p = 1 is  a tangent  plane  to  s}  (8.44) 

It  is  known  that  W is  a real  algebraic  variety  whose  degree  is 
the  class  number  of  S [7,28].  Moreover,  the  relation  of  S and  W 
is  symmetric:  S is  also  the  polar  reciprocal  of  W.  It  is  clear 

that  if 

N(p)  = the  set  of  all  exterior  unit  normals  to  S at  p (8.45) 

then 

W={x=(p*N(p))~1N(p):  p£S)  (8.46) 

Now  the  group  velocity  Vg(p)  = VpA|c(p)  is  normal  to  S at  each 
p € S - Z<*.  Moreover,  p*  VpXjc(p)  = Xj,(p)  = 1 for  such  points  p 
by  (8.39)  and  the  homogeneity  of  Xj^(p)  . Hence 

(x  = vg(p)  H VpXk(p):  p e S - Z^}  C W (8.47) 

for  k = 1,2,3. 


8.10  The  polar  reciprocal  map  T:  S W 

This  is  the  map  defined  in  (8.46);  i.e., 

T(p)  = (p  • N(p) ) 1 N(p)  for  all  p £ S (8.48) 

As  indicated  above,  N(p)  is  not,  in  general,  single  valued.  It 
follows  that  T may  be  neither  single-valued  nor  injective.  How- 
ever, it  was  shown  in  [49]  that  if 

Zg  = set  of  singular  points  of  S 

> (8.49) 

Zy  = set  of  singular  points  of  W j 

Zg  = T-1  Z^,  Z^  = TZg  (8.50) 

ZS  * ZS  U ZS’  ZW  = ZW  U ZW 


(8.51) 
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then  Zg  and  are  sub-varieties  of  dimension  <1  and 

T is  bijective  and  analytic  from  S - Zg  to  W - Zy  (8.52) 


8.11  Examples 

The  equation  (8.37)  for  the  slowness  surface  of  a crystal 
contains  21  independent  parameters  in  the  most  general  case  (tri- 
clinic crystals).  Hence  a great  variety  of  slowness  surfaces  are 
possible.  Crystal  symmetries  may  reduce  the  number  of  parameters. 
The  slowness  surfaces  of  the  various  symmetry  classes  have  been 
studied  by  many  authors.  Thorough  discussions  and  examples  may 
be  found  in  [3]  and  [24]  where  specific  numerical  information  on 
the  stress-strain  tensors  of  real  crystals  may  also  be  found. 

Here  two  examples  will  be  described  briefly  to  show  the  kind  of 
surfaces  that  may  occur. 

Cubic  crystals.  In  this  case  symmetry  reduces  the  number  of 
independent  parameters  to  3 and  the  equation  for  S can  be  written 
[24] 

3 P2 

l J - 1 (8.53) 

j=1  a-b|p|2-cp* 

£in 

Of  course,  the  positive  definitenes  of  c^  imposes  certain  numer- 
ical restrictions  on  a,  b and  c.  Equation  (8.53)  represents  a 
surface  of  degree  6 which  is  irreducible  except  for  special  param- 
eter values. 


Hexagonal  crystals.  In  this  case  symmetry  reduces  the  number 
of  independent  parameters  to  5.  Moreover,  S is  necessarily  a sur- 
face of  revolution  and  reduces  to  two  components  whose  equations 
can  be  written  [24] 


!(pf  + Pz>  + b 2 P 3 


(8.54) 


P1+P2 

c2 -d2 )pj 2+e(p2+p2) 


eL 

-d2|p| 


2+fp2 


= 1 


(8.55) 


(where  a,  b,  c,  d,  e and  f can  be  expressed  in  terms  of  5 indepen- 
dent parameters).  The  two  equations  have  degrees  2 and  4,  respec- 
tively. These  surfaces  of  revolution  can  be  visualized  from  their 
traces  on  the  p1,p3-plane;  see  [24,p.99]  for  a graph  of  such  an  S 
and  the  corresponding  W.  It  is  seen  that  in  the  example  Zg  consists 
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of  2 circles  and  2 points  lying  in  S while  Zy  consists  of  8 circles 
and  2 points  lying  in  W. 


8.12  Asymptotic  wave  functions  for  crystals 

It  was  shown  in  [44]  that  the  equations  (8.1)  for  acoustic 
waves  in  crystals  have  asymptotic  wave  functions  of  the  form 

v(0) 

v (t,x)  = l F(x*s^ (0)  - t,s^ (6))/ |x|  , 

a=l  (8.56) 

x = |x  1 9 

where 

s(a)(6)  e S,  a = 1,2, • •• ,v(0)  (8.57) 

is  the  solution  set  of  the  equation 

N (s)  = 0 (8.58) 

Thus  s(a)  defines  the  multivalued  inverse  of  the  Gauss  map  N of  S. 
The  principal  properties  of  vco(t,x)  are  described  by  the  following 
theorem  whose  proof  is  contained  in  [44]. 


8.13  Theorem 

For  each  h G 3f  there  exists  a unique  F:  R x s •>  Ci3  such  that 

v°°(t,*)  G 3f  for  all  t G R (8.59) 


OO 

t -*■  v (t,*)  G K is  continuous  for  all  t G R 


(8.60) 


Iv  (t,*)llj^  < Cllh'l^  where  C is  independent  of  h and  t (8.61) 

OO 

Finally,  v is  an  asymptotic  wave  function  for  v(t,*) 

= exp  (-itA1/,2)h: 


lim 

t->co 

Moreover, 
in  [49]. 


v ( t , • ) - v ( t , • ) II  = 0 (8.62) 

(a) 

explicit  constructions  of  s (0)  and  F(t,s)  are  given 
In  the  present  case  they  take  the  follovring  form. 


i 


| 

I 


V 
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8.14  Construction  of  s^°^(9) 


The  construction  consists  of  two  steps. 

(0) , a = l,***,v(0)  is  the  intersection  of 

(8.63) 

W - Zy  and  the  ray  from  0 along  0 
s(Ct)(0)  = T_1x(a)(0)  e S - Zs  (8.64) 

fry) 

Note  that  this  defines  s'  (0)  for  all  0 outside  of  Lhe  null  set 


Z°  = {0:  x = 1 x J 0 G Zw)  C S2  = {0:  J 0 J = 1}  (8.65) 


8.15  Construction  of  F(x,s) 

F is  calculated  from  h = v ((},•)  G 3C  by  the  rule 


f °° 

F(x,s)  = (27r)“^2  W(s)  eixX  h(As) AdA 

> 0 

(8.66) 

where 

VCs)  = >(s)  |K(s)  |~1/2  |T(s)!-1  P(s) 

(8.67) 

iKs)  = exp  (i  J ( p~ ( s ) - p+(s))} 

(8.68) 

+ 

p (s)  = the  number  of  principal  curvatures  of 
at  s which  are  ■<  0. 

(8.69) 

K(s)  = Gaussian  curvature  of  S at  s 

(8.70) 

P(s)  = orthogonal  projection  of  C3  onto  the 
eigenspace  for  the  eigenvalue  p = 1 of  A(s) 
(s  6 S) 

(8.71) 

It  is  shown  in  [44]  that  f(s)  is  defined  for  all  s G S - Zg . In 
particular,  the  parabolic  points  of  S lie  in  Zg.  The  integral  fo 
F need  not  converge  pointwise,  but  it  converges  in  the  Hilbert 
space  3f(S)  with  norm  defined  by 
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I |f(t,s) | 2 |k(s)T(s) | dSdT 
•*S 


(8.72) 


Moreover,  the  operator  0:  JC  -*•  JC(S)  defined  by  0h  = F is  an 

isometry. 


8.16  Propagation  in  non-uniform  crystals 

The  method  developed  in  [42]  and  section  4 can  be  applied  to 
local  perturbations  of  uniform  crystals.  Eigenfunction  expansions 
for  non-uniform  crystals,  and  more  general  systems,  have  been 
given  by  G.  Nenciu  [25]. 
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